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An ab-initio calculation scheme for finite nuclei based on self-consistent Green's functions in the 
Gorkov formalism is developed. It aims at describing properties of doubly-magic and semi-magic 
nuclei employing state-of-the-art microscopic nuclear interactions and explicitly treating pairing 
correlations through the breaking of U{1) symmetry associated with particle number conservation. 
The present paper introduces the formalism, necessary to undertake applications at (self-consistent) 
second-order using two-nucleon interactions, in a detailed and self-contained fashion. First applica- 
tions of such a scheme will be reported soon in a forthcoming publication. Future works will extend 
the present scheme to include three-nucleon interactions and implement more advanced truncation 
schemes. 
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I. INTRODUCTION 

Over the last decade the reach of ab-initio nuclear 
structure calculations has extended up to the region of 
medium-mass systems. Despite the significant progress 
from both theoretical and computational points of view, 
methods as coupled-cluster (CC) U, in-medium similar- 
ity renormalization group (IMSRG) [H or Dyson self- 
consistent Green's function 0] (Dyson-SCGF) are how- 
ever currently limited to a few tens of doubly-closed shell 
nuclei. Neighboring nuclei with ±1 or ±2 nucleons can 
also be reached with particle attachment or removal for- 
malisms [3, Q ■ While improving further the convergence 
of such existing schemes, it is essential to extend their 
reach and the intrinsic predictive character of ab-initio 
methods to truly open shell nuclei. One way of doing 
so involves the development of mult i- reference schemes, 
such as e.g. multi-reference CC Q or valence-space shell- 
model based on microscopic inputs from the ab-initio 
calculation of a doubly-closed shell core nucleus of ref- 
erence 0, Q . Alternatively, one may prefer to keep the 
simplicity of a single-reference method. This requires 
however, in any of the approaches mentioned above, to 
formulate the expansion scheme around a vacuum that 
can tackle Cooper pair instabilities, e.g. to build the cor- 
related state starting from a Bogoliubov vacuum that al- 
ready incorporates zeroth-order pairing correlations. The 
objective of the present work is to realize the latter pro- 
gram within the particular frame of SCGF theory. 

As alluded to above, SCGF methods are being suc- 
cessfully applied to the study of nuclear systems. Over 
the last two decades considerable progress has been made 



in the development of suitable formalisms and computa- 
tional algorithms both for finite nuclei and infinite nu- 
clear matter In infinite systems, bulk and single- 
particle properties are typically computed through the 
resummation of particle-particle (pp) and hole-hole (hh) 
ladder diagrams, i.e. in the self-consistent T-matrix 
approximation, that tackles short-range correlations in- 
duced by the hard-core of conventional nucleon-nucleon 
(NN) interactions. Results have been obtained at zero 
and finite temperature for both symmetric and pure neu- 
tron matter based on various conventional NN poten- 
tials [gj-fllj. Recently, microscopic three-nucleon (NNN) 
forces have been incorporated [l2|, [T^] . There have been 
also attempts to take into account nucleonic superfiuidity 
through the consistent treatment of anomalous propaga- 
tors 

In finite systems the most advanced SCGF calcula- 
tions feature the Faddecv random-phase approximation 
(FRPA) technique, which allows the simultaneous inclu- 
sion of pp, hh and ph excitations, together with inter- 
ferences among them [l^ [l3 • By employing a G- matrix 
resummation of scattering diagrams not included in the 
chosen model space it is also possible to use interactions 
with strong repulsive cores |18| . An important charac- 
teristic of the FRPA expansion is that it is based on 
combining one- and many-body propagators, each one 
representing different experimental processes including 
nuclear excitations and transfer of one or two nucleons. 
The method has therefore been applied to a variety of 
jroblems including the quenching of spectroscopic factors 
, anharmonic excitations |20|, two-nucleon knockout 
[22 1, and the derivation of optical potentials [23lf23. 
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At the moment, applications can access all doubly-closed 
shell nuclei up to ^^Ni together with neighboring systems 
with ±1 or ±2 nucleons @, [11 . 

In the present work SCGF calculations of finite nu- 
clei are implemented within the Gorkov scheme, allow- 
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ing for an explicit treatment of nucleonic superfluidity. 
Suitable numerical techniques are developed in order to 
perform systematic studies of doubly-magic and semi- 
magic medium-mass nuclei as will be soon reported on in 
a forthcoming publication [2^, referred to thereafter as 
Paper II. 

One of our goals is to be able to tackle various types of 
nuclear interactions, in particular chiral potentials based 
on effective field theory (EFT) [2^ and low-momentum 
potentials obtained through the further ap plication of 
renormalization group (RG) techniques [27[. There arc 
also yet unanswered fundamental questions as to what 
microscopic processes are res pon sible for the superfluid 
character of open-shell nuclei |28l - [37| . This is one among 
several long-term objectives of the project to provide a 
fully ab-initio answer to such questions. The present 
work eventually relates as well to the long-term devel- 
opment of so-called non- empirical energy density func- 
tionals (EDFs) HH, [H, [H . There exist on-going efforts 
to construct nuclear EDFs starting from underlying nu- 
clear interactions, with the main goal of improving the 
predictive power away from known data that is rather 
poor for existing phenomenological EDF parameteriza- 
tions. The connection with NN and NNN interactions 
is typically obtained by means of density matrix expan- 
sion (DME) techniques [13 and many-body perturbation 
theory, which allow for the construction of schemes that 
can be systematically tested and improved order by order 
in the interaction [40-H]. In this regard, recent develop- 
ments and applications of low-momentum potentials |45l - 
l47j , which exhibit a more perturbative nature than tradi- 
tional nuclear interactions, are instrumental. Schemes to- 
wards non-empirical EDFs, however, are presently avail- 
able only in their first stages. Their development ne- 
cessitates a comparison with fully microscopic methods 
that can provide useful benchmarks over which EDFs pa- 
rameterizations can be tested and improved. In this con- 
text, SCGF techniques represent a valid ab-initio method 
of reference. In particular, as the breaking and restora- 
tion of symmetries (e.g. translational, rotational, particle 
number, ...) is central to nuclear EDF methods, it is cru- 
cial to develop an approach that includes and exploits the 
same concept [48|, which is the case of Gorkov-Green's 



function method regarding particle-number symmetry. 

The present paper aims at providing a detailed account 
of Gorkov's formalism that is eventually applied in Pa- 
per II. Given that the present work is the first nuclear 
structure application of ab-initio Gorkov self-consistent 
Green's function method, it is relevant to provide a self- 
contained account of the formalism expressed in a dis- 
crete basis, which is suited to finite nuclear systems. The 
present formalism is further specified to second-order in 
the self-energy expansion and formulated in terms of NN 
interactions only. The extension to more advanced trun- 
cation schemes and to the inclusion of NNN forces is post- 
poned to future works. 

The paper is organized as follows. Scctionllllintroduces 
the general form of the nuclear Hamiltonian employed in 



the present work while Sec. IIIII defines generic features of 
Gorkov's formalism. Section |IV] discusses Gorkov's equa- 
tion of motion under the form of an energy-dependent 
eigenvalue problem before computing normal and anoma- 
lous self-energies at second order and rewriting Gorkov's 
equation under the form of a more convenient energy- 
independent eigenvalue problem. Further details regard- 
ing the extraction of observable arc provided in Sec. |Vl 
while Sec. |VT] discusses the conserving character of the 
employed truncation scheme. Finally, conclusions are 
given in Sec. IVIII followed by several appendices com- 
plementing the body of the paper with relevant technical 
details. 



II. NUCLEAR HAMILTONIAN 

A. Single-particle basis 

Let us first introduce the particular labeling of single- 
particle states that will be used throughout the work. 
We consider a basis {aj,} of the one-body Hilbert space 
Hi that can be divided into two blocks according to the 
value (or more precisely to the sign) of an appropriate 
symmetry quantum number. To any state a belonging to 
the first block, one can associate a single-particle state 
a belonging to the second block and having the same 
quantum numbers as a, except for the one differentiating 
the two blocks. Typically there exists an anti-unitary 
transformation T, leaving the Hamiltonian invariant and 
connecting, up to a phase r]a, state a with state a. With 
that in mind one defines a basis {aj^}, partner of the 
initial one {aj,}, through 



al {t) = riaal {t) , aa (t) = rjaCid {t) 



(1) 



which corresponds to exchanging the state a by its part- 
ner d up to the phase r]a- By convention d = a with 

As discussed in Sec. O T will eventually be specified 
as the time-reversal transformation in our applications to 
even-even nuclei with = 0^ ground states. 



B. Hamiltonian 

Let us consider a finite system of N fermions interact- 
ing via a two-body potential 1/^^. The corresponding 
Hamiltonian can be written as 



(2!)2 ^"'""^ alaladUc (2) 



E 

ibcd 



^ In the following, a tilde (~) always refer to the sole quantum 
numbers of the opposite block while barred (~) quantities involve 
an additional phase factor r^^^g. 
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where 



or as a straight two-body operator 



= {a\T\b) 



(3) 



is the matrix element of the kinetic energy operator T 
and 

Vabcd = (a&|V^^^|cd) (4) 
= (l:a; 2:b\V^^ \l:c; 2:d) - (l:a; 2:b\V^^ \l:d; 2:c) 

is the antisymmetrized matrix element of v'^'^ expanded 
in terms of direct-product states, denoted by |1,2). 
Whenever a superscript or subscript index x appears in 
a matrix element (jS)) or Q, the associated annihilation 
(creation) operator is to be intended of the form ([T]), i.e. 
cix (aj.). It follows that, e.g.. 



Tab = VaVb ia\T\b) 



(5a) 



Vabcd = VaVc{ib\V''''\cd) 



(5b) 

etc. As discussed in Appendix IC 1[ properties of operator 
T leads to the following useful relations 



T*, 



b ' 



(6a) 



T. 



(b) _ 



N 



2M 



(11) 



Here pi represents the momentum of the z-th nucleon, M 
the nucleon mass and N is the particle number operator. 
In theories that do not conserve particle number, N can- 
not be replaced by its eigenvalue such that expressions 
and ([TT|) are not equivalent. Considering a series 
expansion in N~^, it could be shown [i^ that form PH)) 
displays the correct power counting and should be there- 
fore employed in calculations over Fock space. 

In the following, we consider Hamiltonian (j9|) with 
choice ([TU]) at first order in N^^, i.e. N^^ is replaced 
by its average value iV~^. For simplicity, and unless oth- 
erwise stated, we denote in the following H,-ci by H such 
that T actually embodies the one-body part of T^^^ (first 
term in Eq. (fTU]) ) and such that incorporates the 

two-body part of T^^^ (second term in Eq. ([TU])). 



III. GORKOV FORMALISM 



^abcd 



abed 



(6b) 



A. Standard propagator and superfluid systems 



C. Centre-of-mass correction 

In the study of a N-body self-bound system, a sep- 
aration can be made between the motion of its centre- 
of-mass and the motion of the nucleons relative to it. 
Specifically, the N-body Hamiltonian ([2]) can be divided 
into 



(7) 



where i?cm represents the centre-of-mass kinetic energy 
and the internal Hamiltonian iJrei does not depend on 
centre-of-mass coordinates. Eigenfunctions of fftot can 
therefore be expressed as products of eigenfunctions of 
Hem and eigenfunctions of -ffroi- Consequently, the en- 
ergy is the sum of centre-mass energy Ecm and internal 
energy E^^i 



Et, 



Er: 



E, 



rcl 



(8) 



Nuclei being self-bound objects, one is interested in the 
translationally invariant, internal Hamiltonian H-^d and 
the corresponding energy E^^i- Subtracting the (known) 
centre-of-mass kinetic contribution from the total Hamil- 
tonian, one indeed works with the internal Hamiltonian 



^^rcl — ^^tot — Hem — Tic 



V 



NN 



(9) 



The internal kinetic energy can be expressed either as a 
sum of one- and a two-body operators 



T. 



(a) 



rcl 



1 



1 

N 



y- 



2M N 



M 



(10) 



of 



Let us consider the N-body ground-state |^'q ) solution 



H 



(12) 



with the lowest eigenvalue Eq . The fundamental object 
of Green's function theory is the one-body propagator 
defined as 



\t,t')^{^^\T{aa{t)al{t')}\^SJ^), (13) 



where the operator T orders a and according to their 
time argument (larger times to the left) and where an- 
nihilation and creation operators are in the Heisenberg 
representation 

Uait) = ai^\t) = exTp[iHt\ aaexp[—iHt\ , (14a) 



4(t) 



AH) 



(t) = cxp[iHt] 4 cxp[-iHt] . (14b) 



The knowledge of Q enables the computation of expecta- 
tion values of all one-body operators plus the two-body 
ground-state energy, i.e. the expectation value of the 
Hamiltonian if only two-body forces are considered. One 
can define two-, three-, A^-body propagators in a simi- 
lar way, in order to evaluate up to A-nucleon observables. 

Green's functions' equations of motion take the form 
of a set of A^ coupled integro-differential equations, each 
of them involving (i — 1)-, i- and (i-l-l)-body propaga- 
tors. In order to compute the one-body propagator, one 
can as well derive a perturbative expansion that trans- 
lates into an infinite series of diagrams. Both approaches 
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provide systematic ways of approximating the exact solu- 
tion. The connection between the diagrammatic expan- 
sion and the equation of motion for Q leads to the defini- 
tion of the (irreducible) self-energy E and the derivation 
of Dyson's equation 

Gir' = sir' + E ^i^"^ ^^r"^ ' (15) 

cd 

where t/*-*^-* is the one-body propagator of the unperturbed 
system associated with a one-body Hamiltonian Hq of 
choice. 

The scheme is in principle exact, i.e. if one can com- 
pute the perturbative expansion up to infinite order. Ap- 
proximations are introduced by including only a certain 
class or subset of terms in the computation of the self- 
energy. Such a subset is chosen according to a hierarchy 
between the various types of diagrams whose rationale 
depends on the system under consideration. The valid- 
ity of the standard perturbative expansion, however, is 
not always guaranteed. In particular, nuclear interac- 
tions inducing strong pairing correlations between con- 
stituents of the many-body system make the usual ex- 
pansion inappropriate for the large majority of nuclei. 
The breakdown of the perturbative expansion is signaled 
by the appearance of (Cooper) instabilities, which occur 
when summing up certain classes of diagrams and point 
to the necessity of developing an alternative diagram- 
matic method |50l - [53 |. 

B. Auxiliary many-body problem 

In the presence of Cooper instabilities, one can de- 
velop an alternative expansion method accounting in a 
controlled fashion for the appearance and destruction of 
condensed nucleonic pairs. 

Instead of targeting the actual ground state j'l'^ ) of the 
system, one considers a symmetry breaking state |4'o), 
i.e. a wave packet, defined as a superposition of actual 
ground states of (A^— 2)-, A^-, (A^-|-2)-, ... particle systems, 
i.e. 

even 

|*o)^E^^lO' (16) 

N 

where cn denote unknown complex coefficients. The sum 
over even particle numbers is said to respect the (even) 
number-parity quantum number. Together with such 
a state, one considers the grand-canonical-like potential 
VL = H — fiN, with fi the chemical potential, in place of 
hB The state l^-o) is chosen to minimize 

no = i^oimo) (17) 



^ Let us remark that the analogy with a grand-canonical ensemble 
holds only on a formal level, as here eigenstates of Q are pure 
states IV'q'), \4>q^'^), ■■■ while the admixture or symmetry break- 
ing state l^o) is a pure state as well. In other words, one is not 
introducing a statistical density operator to describe the system. 



under the constraint 

AT = (*o|A|*o) , (18) 

i.e. it is not an eigenstate of the particle number op- 
erator but it has a fixed number of particles on average. 
Equation (|17p , together with the normalization condition 

even 

(*o|*o) = El^^l' = l' (19) 

N 

determines the set of coefficients cn ■, while Eq. (fT8|) fixes 
the chemical potential /i. 

By targeting |^'o), the initial problem that aimed at 
describing the many-body system with A^ nucleons is 
replaced with an auxiliary problem, whose solution ap- 
proximates the initial one. The validity of such an ap- 
proximation resides in the degeneracy characterizing the 
ground state of the system. The presence of a condensate 
(ideally) implies that pairs of nucleons can be added or 
removed from the ground-state of the system with the 
same energy cost, independently of A^. Such an hypoth- 
esis translates into the fact that the binding energies of 
the systems with A^, N±2, A^±4, ... particles differ by 2//; 
i.e. the idealized situation considered here corresponds 
to the ansatz that all ground states obtained from the 
system with A" nucleons by removing or adding pairs of 
particles are degenerate eigenstates of such that their 
binding energies fulfill 

... « - < « < - «... « 2m , (20) 

with M independent of A. If the assumption is valid, 
the energy obtained by solving the auxiliary many-body 
problem provides the energy of the initial problem as 

l^o = El^^'l'^^' ^^0^-*"^' (21) 

N' 

which follows from Eqs. (HH) and (^0)) . 

C. Gorkov propagators 

In order to access all one-body information contained 
in 1 5*0), one must generalize the one-body propagator 
defined in Eq. by introducing additional objects that 
account for the formation and destruction of pairs. One 
thus defines a set of four Green's functions, known as 
Gorkov propagators [s^l , through 

I G]^,{t, t') ^ (*o|T [aa{t)al{t')] \^o) , (22a) 
I Gllit, t') = (*o|T {aa{t)ah{t')} |*o) , (22b) 
zGllit,t')^{^o\T{aiit)al{t')}\^o), (22c) 
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^Gll{t,t')^{^o\T{ai{t)a,{t')}\^o}, 



(22d) 



where single-particle operators associated with the part- 
ner basis are as defined in Eq. ((T]) and where the modified 
Heisenberg representation is introduced through 

aa{t) ~ a^a^\t) = cxp[iflt] Qa cxp[—iflt] , (23a) 

4W = aa^'W ^ ^cxp[int]alexp[-int] . (23b) 

Besides the time dependence and quantum numbers a 
and b identifying single-particle states, Gorkov propaga- 
tors j,^^ carry two labels gi and 52 that span Gorkov's 
space. When gi = 1 {gi = 2) a particle is annihilated in 
the block of a (created in the block of 5) and vice versa 
for (72; i-e. 32 = 1 (.92 = 2) corresponds to a second par- 
ticle created in the block of b (annihilated in the block 
of b). Green's functions G^^ and G^^ are called normal 
propagators while off-diagonal ones, G^^ and G^^, are 
denoted as anomalous propagators. 



E. Energy representation 

For most applications it is convenient to transform 
the propagators from time to energy representation. In 
systems at equilibrium governed by a time-independent 
Hamiltonian, one-body Green's functions depend only 
on the difference of their two time arguments, i.e. 
Gab{t,t') = Gab{t — t'). Gorkov propagators in the en- 
ergy domain arc thus obtained through the Fourier trans- 
formation 



Gob(w) 



d{t-t')e''^^'-''^Gab{t~t') ■ (28) 



The energy representation is more suitable to analyzing 
the physical content of singlc-particlc propagators, as will 
become clear in the following. 



F. Gorkov's equations 



D. Nambu's matrix formalism 

Gorkov's propagators can be conveniently grouped into 
a matrix representation, introduced by Nambu |54l | . First 
one defines the two-component vector 



A„(f) 



aa{t) 

4(0 



and its self adjoint 

A\{t)^{a\{t) -aa{t)) , 



(24a) 



(24b) 



denoting generalized annihilation and creation operators. 
Their components fulfill the anti-commutation relations 

[Al^{t),Af\i)] = 5,,,Jab. (25) 

where the extra label labels the rows (columns) of the 
annihilation (creation) vector operator. One can then 
write the four propagators (|22p in the matrix form 



(26) 



iGabit.t') = («'o|T{Aa(OAt(i')} l^-o) 
^G]i{t,t') Gll{t,t'y 

Gllit,t') Gll{t,t') 



where the time ordering operator acts separately on each 
element of the Gorkov's matrix AA^. In general, any 
object 0^1^^ defined in Gorkov's space can be put into 
such a matrix form 



Oabit,t') = 



'oll{t,t') oil{t,t')' 
ofbit.t') oll{t,t' 



(27) 



with gi and (72 labeling respectively the rows and the 
columns of the matrix. 



In the standard case, the derivation of the equations of 
motion and the formulation of a diagrammatic expansion 
for the one-body propagator lead to defining the irre- 
ducible self-energy and Dyson's equation, through which 
the propagator of the interacting system can actually be 
computed. One proceeds similarly in the Gorkov formal- 
ism. The first step consists in separating the Hamiltonian 
into an "unperturbed" one-body part and an interact- 
ing part. This is conveniently achieved by introducing 
an auxiliary, one-body Hermitian potential U taking the 
general form 



U = 



^Uab a\ab - Uab aaal + Uab alal + [7]^, aaab 

ab 



and by defining 

n=.T + U-^iN + v^^ - u 



(29) 
(30) 



The one-body grand potential flu defines unperturbed 
Gorkov propagators G'-^-' in energy representation 
through 



G(")(a;)j =uj-ni 



(31) 



The choice of G'^"'' corresponds to selecting an appropri- 
ate unperturbed ground state which acts as a reference 
vacuum for the application of Wick's theorem, and is 
crucial for the convergence of the perturbative series. In 
particular, one cannot expand the interacting superfiuid 
ground-state |^'o) around a non-superfluid unperturbed 
state, i.e. unperturbed propagators must already contain 
the basic features characterizing the interacting ones. 

The requirement that the unperturbed ground state is 
superfiuid translates into a choice of a flu that breaks 
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particle number, as is evident from the form of the aux- 
ihary potential (P^ . Applying Wick's theorem in the 
derivation of the perturbative expansion, anomalous con- 
tractions naturally appear, and are afterwards identified 
with anomalous Gorkov propagators. 

Once the unperturbed ground state is defined, one 
writes down the perturbative series for the interact- 
ing propagator G and defines normal and anomalous 
one-line irreducible self- energies. Self-consistency is ob- 
tained by computing self-energy diagrams in terms of 
fully dressed propagators G and by only retaining skele- 
ton diagrams, i.e. diagrams with no self-energy insertions 
(see App. lB l]) . Working in the energy representation, the 
four irreducible self-energies read 



Sab(a;) 



(32) 



and can be divided into a proper part and a contribution 
coming from the auxiliary potential, i.e. 



Sa6(w) = Sa6(a;) - Uq, 



b ■ 



(33) 



Eventually, standard Dyson's equation is generalized as 
set of coupled equations involving the two types of prop- 
agators and self-energies. These are known as Gorkov's 
equations [s^ and read, in Nambu's notation, 



equations once the propagators are dressed with the cor- 
responding self-energies. From a practical point of view 
it is useful to track where the auxiliary potential enters 
and how its cancelation is eventually worked out. This 
point is addressed in Section IIV A[ where the solution of 
Gorkov's equations is discussed. In particular, and since 
such a solution is to be found through an iterative proce- 
dure, one is still interested in choosing a good auxiliary 
potential as a starting point. 

Let us further remark that, as the auxiliary potential 
(Uni) has a one-body character, i.e. it acts as a mean field, 
the search for the ground state of fie/ corresponds to solv- 
ing a Bogoliubov-like problem, as becomes evident when 
writing the unperturbed grand potential in its Nambu's 
form 



[n 



U\ab 



Tab + Uab — M Sab 



U. 



Uab 

'Tab — Uab 



fJ.Sab 



(35) 

In fact a convenient choice for riu is constituted by 
^HFB, i-e. one first solves the Hartree-Fock-Bogoliubov 
(HFB) problem and then uses the resulting propagators 
G^ff^ as the unperturbed ones. Notice that the self- 
energy corresponding to this solution, S^^^ ^ eventually 
differs from the first-order self-energy S^^^ as soon as 
higher orders arc included in the calculation because of 
the associated self-consistent dressing of the one-body 
propagators. 



GabiLo) = Gi",)(c.) +^ Gi°)(L.) i:,rf(w) G,b(a 



(34) 



As Dyson's equation in the standard case, Gorkov's equa- 
tions represent an expansion of interacting or dressed 
one-body normal and anomalous Green's functions in 
terms of unperturbed ones. If the method is self- 
consistent, the final result does not depend on the choice 
of the auxiliary potential, which disappears from the 

I 



G. Lehmann representation 

Let us consider a complete set of normalized eigen- 
states of with no definite particle number 



(36) 



and spanning the Fock space J^. Inserting the corre- 
sponding completeness relation, G^^{t,t') becomes 



Gll{t,t') = -^e{t - f) Y,{^^\aa\^k){^k\al\^o) > 



,i[no-^k]{t-t') 



+ ie{t' - t) ^(*o|aIl*fc)(*fc|aa|*o) ' 



-i[no-nfc](t-t') 



r 



Using the integral representation of the theta function Lehmann representations 
and reading out the Fourier transform, one obtains the 
propagator in energy representation under the form 



{^o\aa\^k){^k\al\^o) 

cj - [rifc - r^o] -I- IT] 

(^o|4l^fc)(^fc|a„|^o) 
w -f [rife - r^o] - irj 



(37) 



GJ^H = E^ '''''' 



= E 



yfe* yk 
a '^b 



w - Wfe + if] oj + uik — ir] 



uj — ujk + iif w + Wfe — z?7 



(38a) 



(38b) 



One can proceed similarly for the other three Gorkov- G'^\{uj) = 

Green's functions and obtain the following set of 



VaK* 



cj — UJk + ir] UJ + uJk — irj . 



(38c) 
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yfc yfc* 
a b 



uj - bJk + iri uj + uJk — iri 



(38d) 



with Gorkov's spectroscopic amplitudes defined as 

U^* ^{^k\al\^o), (39a) 

Vt = {^k\aa\^o), (39b) 

and 

U^*^{^k\al\^o), (40a) 

Vt = {^k\aa\^o), (40b) 
from wlriclr follows thal|f| 

Ui = +i^aUl, (41a) 

V^a=-VaV^a- (41b) 

The poles of the propagators^ are given by tok = fi^ — fio. 
The relation of such poles to separation energies between 
the N-body ground state and eigenstates of the A^±l sys- 
tems is polluted by the breaking of particle number sym- 
metry and is less transparent than for standard Dyson- 
Green's function. Still, the structure of the propagator 
naturally suggests to approximate one-nucleon separa- 
tion energies as 

-/x[(*fe|iV|*fc)-(^ + l)] , (42a) 

+^[(*fc|7V|*fc)-(iV-l)] , (42b) 

where the error associated with the difference between 
the average number of particles in state and the 

targeted particle number ± 1 is taken care of by the 
last term in Eqs. (|42ap and (|42b|) . 

It is useful to introduce a Nambu representation for 
the Lehmann form of the propagators by defining the 
two-component vectors 

X^-t^(vl/fe|At |vl/o) = (Z^f ), (43a) 



(*fc|AJ*o) 



(43b) 



where A and A^ have been introduced in Eq. (p4|) . and 
by writing 



(44) 



Note that vectors (|33]) contain equivalent physics infor- 
mation and are transformed into each other by 



x!: = 



-1 

1 



(45) 



^ Similarly to Eq. ((S}, wc may equivalcntly write Eqs. l|4H l as Uj^ = 
and = -V|. 

* As discussed later on, eigensolutions of Gorkov's equations come 
in pairs (cj^,,— oj^.) such that one should only sum on positive 
solutions in Eq. 13811 . 



H. Symmetry properties 

The four Gorkov propagators and self-energies are not 
independent from each other and can be related through 
certain symmetry operations. Starting from the defini- 
tion of Gorkov Green's functions ((22)) and their Fourier 
transforms ()28|) . one can first prove that 



Glliio) - - rjaVb Glli~u) = - Ggi~u) , (46a) 
Glliio) ^+rjar^bGll{-io) ^+GII{-l,) , (46b) 
Glliio) ^+rjaVbGll{~co) ^+Gll{-co) . (46c) 



Result ()46p is easily derived from Lehmann representa- 
tion (p8|) . together with properties (|4T|) . as 



Glli^) 



yfc yfc* 
^ a ^ b 



El 

[UJ - UJk + IT] UJ +UJk- IT] 



E 



—UJ + ujk — irj —UJ — ujk + ii] 



(47a) 



and 



Gfbi^) = 



^ \ uj — uJk + irj UJ + ujk — iri 



E 



-UJ + UJk — irj —UJ — UJk + irj 



G (-) 



(47b) 



By separating the real and imaginary parts of the poles in 
Eq. (|44p the Gorkov propagator splits into its hermitian 
and antihermitian components 



Cliuj) Aabiuj) 



+ i 



Bab{uj) Dab{uj) 
Dliuj) Bab{uj) 



(48) 

where Aab{uj) and Bab{uj) are hermitian matrices in 
the one-body Hilbert space "Hi. Note that because of 
the presence of an antihermitian component G]^(uj) ^ 
\G^{uj)\* and Gf,{uj) + {Gf,{uj)X- From S and 
(|46E| it follows that 

Kh{^) = -Ai-^{-uj) , 
Bab{^) = -Bi^{-uj) , 

Cabiuj) = +C-^a{-Uj) , 
Dabiuj) - +D-,^{-Lj) . 

Similar symmetry properties are valid for normal and 
anomalous self-energies. Starting from Gorkov's equation 
(|34| and making use of relations (|46)) one can prove that 
the equivalence between Gorkov's equation and its 
conjugate (|A9|) requires 

J:ll{uj) = - i^aVb mi-^) = ~ mi-^) , (49a) 
J:i:i{uj) = +r^aVbm{-^) = +^lli-co) , (49b) 
J:IUuj) = +r,aVbmi-uj) = +m{-u) . (49c) 
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Such properties are general and should be required from 
any truncation scheme used to compute self-energies. At 
first order, they arc confirmed by the explicit evaluation 
of normal and anomalous diagrams in Eq. (j65p . At sec- 
ond order one can check that they are indeed fulfilled by 
expressions fTT]) and (f79| . 



I. Spectroscopic content of Gorkov propagators 

Let us now discuss quantities that arc useful to analyze 
the spectroscopic content of Gorkov propagators. First, 
one defines generalized spectroscopic factors through the 
2x2 Nambu matrix 

F,^^(vI'o|AJvi/,)(vi/,|At|vl/o) (50) 

a 
a 

which is independent of the one-body basis used and 
whose normal components generalize traditional spectro- 
scopic factors for addition and removal of a nucleon 



TZ = 



a 
a 
a 



(51a) 



(51b) 



t|2 



As states l^'o) and j^'fc) do not carry a definite parti- 
cle number, such spectroscopic factors do not possess the 
sharp physical interpretation of the usual ones. Still, and 
although {J^^ ) contains contributions from the addi- 
tion (removal) of a nucleon to (from) systems character- 
ized by different particle numbers, the dominating con- 
tribution remains associated with the addition (removal) 
to (from) the actual targeted ground-state l^'o')- 

Next is Gorkov's one-nucleon spectral function S(cj) 
summing one-nucleon addition S^(cj) and removal S~ (uj) 
components. Such spectral functions are not only 
(energy-dependent) 2x2 matrices in Nambu space but 
also matrices on the one-body Hilbert space Hi- They 
are extracted from the imaginary part of Gorkov's prop- 
agators through 



S+,(^) 



1 



ImGab(tj) 



E '^(^ - ^fc) for w > , (52a) 



where only Wfc = ±(i?^ — m) ^ contribute to the sum. 
Just as for spectroscopic factors, the normal components 
of Gorkov's spectral functions, e.g., 

S'abiE) = S+'\E-f,) ^ E^a '5(i? - E+) , (53a) 



S'^biE) = S;,'\E-t,) = E ^a* 5{E - E-) , (53b) 

k 

generalize standard particle and hole spectral functions. 
The normal one-body density matrix can be extracted 
by integrating the normal part of the removal spectral 
function '5'~{,^^, 

Pab = (^olabflal^o) (54a) 
du:S~,^\oj) (54b) 

= Y.^bVa\ 

k 

whereas the anomalous density matrix is obtained as 

(54c) 

du:S-,'-\Lo) (54d) 

k 

Information contained in the spectral function can be 
characterized by computing its various moments, i.e. 

/oo 
dw^"S,fc(c^), (55) 
-oo 

which are (energy-independent) matrices in Nambu space 
and on "Hi . Making use of anti-commutation relation ([25]) 
one can derive a sum rule for the generalized spectral 
function that directly relates to its zcroth momenl[£|, i.e. 



Pab = (*o|af)aa|*o) 

12/ 



M. 



(0) fflS2 



{^o\{Ai\Al^^}\^o)^5,,gJ,b. (56) 



The usual sum rule associated with the normal part of 
the spectral function is recovered from Eq. ([56| as 



-foo 



di.[S-,'\i.) + S+'\uj)]=6ab.. 



(57) 



showing that diagonal matrix elements S^l^{iu) are noth- 
ing but probability distribution functions associated with 
the probability to remove/add a nucleon from/to the 
ground state from/on a given single-particle state a and 
leave the residual system with a missing energy w. Equa- 
tion ([57)) simply states that such a probability integrate 



s:,(^) 



■ ImGat(w) 



EY^Y,^^^(^ + '^fe) forc^<0, (52b) 



^ As discussed in Sec. IV Bl the first moment M^^^ of Gorkov spec- 
tral function gives access to effective single-particle energies |55|| . 



9 



to 1 when scanning missing energies from — cx) to +00. A 
new sum rule associated with anomalous spectral func- 
tions can also be deduced from Eq. (|56|) as 



+00 



du:[S~,'\u;) + S+'\i,)]^0. (58) 



IV. GORKOV'S EQUATIONS 

We now proceed to a form of Gorkov's equations al- 
lowing for a direct numerical implementation. 



A. Energy-dependent eigenvalue problem 



Last but not least, one introduces the spectral strength 
distribution (SSD) through Sp(£') = Trn, [S{E-fi)], 
which reads as 



Sp{E) - 5] F+ S{E - E+) + F,T 5{E - E^) . (59) 



Let us first transform Eq. ()34p into an eigenvalue 
equation for amplitudes and V'^, along with a nor- 
malization condition for those amplitudes. Multiplying 
Gorkov's equation (p4|) by (w — uok), the pole at cj = +LOk 
is extracted by taking the limit w — > w^, such that substi- 
tuting Lehmann representation (|44p for G and operator 
form ([3T|l for G'"^, one obtains 



The SSD is a 2x2 matrix of energy-dependent functions 
and is independent of the single-particle basis used to 
compute it. Its normal part Sp{E) = Sp^^(£') reads as 



Sp{E) = - Et) + .F- 8{E -£;-), (60) 



cd 



Multiplying both sides by (w — ^v^ea summing over 
a yields 



and provides the probability to leave the system with 
relative energy i? by adding/removing a nucleon to/from such that Eqs. (|33p and (|35|1 finally allows writing the 



ground state j^fc 



matrix eigenvalue equation 



E 



S21(C.) -r„b + M<5ab + S22(c.) 



(61) 



whose solutions are amplitudes {TU^ ,\>^) and associated pole energy uiu- Equivalently, computing the residue at 
w = — Wfe leads to 



E 



E21(,,) -r,, + /i^,fc + Ef,(c.) 



_a 



(62) 



Notice that the latter relationship can be also obtained 
from the conjugate of Eq. ()6ip by using properties of 
Gorkov amplitudes and self-energies. Equations (|6ip 
or (|62p and their solutions are independent of auxiliary 
potential [/, which canceled out. This leaves proper 
self-energy contributions only, which eventually act as 
energy-dependent potentials. The self-energies depend 
in turn on amplitudes and such that Eqs. (pT|) 
or ((62|) must be solved iteratively. At each iteration the 
chemical potential must be fixed such that Eq. is 
fulfilled, which translates into the necessity for amplitude 
V to satisfy 



(63) 



where pah is the (normal) one-body density matrix (|54ap . 



As demonstrated in App. \^ the spectroscopic ampli- 
tudes solution of Eq. (|6ip or fulfill normalization 
conditions 



Eixsr=i+EXa^ 



l9Sa&('^) 



EiY^r = i+EYr 



9w 

9Sab(w) 



(9w 



, (64a) 



Y,^^ , (64b) 



where only the proper self-energy appears because of the 
energy independence of the auxiliary potential. 
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B. First-order self-energies 

In Fig. [U first-order diagrams contributing to normal 
and anomalous self-energies are displayed. Diagrammatic 
rules appropriate to the computation of Gorkov's prop- 
agators and for the evaluation of self-energy diagrams 
are discussed in App. [B] while the $-derivability of the 
presently used truncation scheme is addressed in Sec. lVIl 

The four first-order self-energies diagrams are com- 
puted in Eqs. ((BT8| . ((B20| . (|B2^ and ((B24| . and read 





— + ^ Vacbd Pdc — 

cd 

— — V- n* ~ 

/ J *bddc Pcd 

cd 




+Al, (65a) 


„22(1) 
^ab 


-A*- 
^^db ' 


(65b) 


Vl2(l) 
^ab 


~ 2 $Z ^abcd Pcd = 
cd 




(65c) 


V21(l) 
^ab 


~ 2 ^ bacdPcd — 




(65d) 



where the normal (pa&) and anomalous [fiab] density ma- 
trices have been defined in Eqs. (|5T 



• 





FIG. 1. First-order normal E^^ "-^^ (left) and anomalous 
j]2i (1) (■j.jgjj^^ self-energy diagrams. Double lines denote self- 
consistent normal (two arrows in the same direction) and 
anomalous (two arrows in opposite directions) propagators 
while dashed lines embody antisymmetrized matrix elements 
of the NN interaction. 



C. HFB limit 

Neglecting higher-order contributions to the self- 
energy, Eqs. (|5T|) and ([55)1 combine to give 



E 

b 



T, 



ah 



Kb 



pSab 



'l-ab 
ah ah 



P Sab 



(66) 



which is nothing but the HFB eigenvalue problem in the 
case where time-reversal invariance is not assumed. In 
such a limit, U'^ and V'^ define the unitary Bogoliubov 
transformation (56j according to 

aa = Y.U^:i3k + VtPi. (67a) 
fe 

ai^J2^tPi + V^f3k. (67b) 

fe 

Moreover, normalization condition (j64b|) reduces in this 
case to the well-known HFB identity 

Y.K\' = Y.H\'+Y.K\' = i- (68) 

a a a 

Let us now stress that, despite the energy independence 
of first-order self energies, some fragmentation of the 
single-particle strength is already accounted for at the 
HFB level such that one deals with quasi-particle de- 
grees of freedom. In particular one can deduce from Eq. 
([68|) that (generalized) spectroscopic factors defined in 
Eq. ([5T|) are already smaller than one. Such a frag- 
mentation is an established consequence of static pairing 
correlations that are explicitly treated at the HFB level 
through particle number symmetry breaking. 

Finally, let us underline again that, whenever higher 
orders are to be included in the calculation, first-order 
self-energies (j65p are self-consistcntly modified (in par- 



ticular through the further fragmentation of the quasi- 
particle strength) such that they do not correspond any- 
more to standard Hartree-Fock and Bogoliubov poten- 
tials, in spite of their energy independence. They ac- 
tually correspond to the energy-independent part of the 
(dynamically) correlated self-energy. 



D. Second-order self-energies 

Let us now discuss second-order contributions to nor- 
mal and anomalous (irreducible) self-energies. In Figs. 




FIG. 2. Second-order normal self-energies E^^ ^-^ ' (left) and 
{2 ) (right). See Fig. [T] for conventions. 

[5] and [3] the four types of normal and anomalous self- 
energies are depicted. The evaluation of all second-order 
diagrams is performed in App. [B] Before addressing 
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and 



FIG. 3. Second-order anomalous self-energies E^^ ' (left) 
and (2 ) (right). See Fig. [l]for conventions. 



their expressions, let us introduce useful quantities 



and 



'pkik2k3 
' a 

-T-)fclfc2fc3 

\fkik2k'i 

Skik2k'i 
a 

Qkik2k3 



ijk 



1^3 ^2 



(69a) 
(69b) 



ijk 
ijk 



ijk 



1^3 ^2 



(70a) 
(70b) 



ijk 

Y.^akijl^t^-'Vt =Nt^^-^- , (70c) 

ijk 



in terms of which second-order self-energies are expressed 
below. Using relations (|^T|) one shows that 



\jklk2k3 




(71a) 


■pkik2k3 
' a 


— n ■T>^i^2k3 

— 'la 1 a ' 


(71b) 


'Ukik2k3 


„ ■7pfclfc2fc3 

— 'la l^a ' 


(71c) 



and 



\rkik2k3 _ 


Kfkik2k3 


(72a) 


Qkik2k3 _ 


r)kik2k3 
'la ' 


(72b) 


^kik2k3 


„ Cfclfc2fc3 


(72c) 



Given that V and TZ can be obtained from Ai through 
odd permutations of indices {A;i,/c2,fc3} and taking into 
account the symmetries of interaction matrix elements, 
one can prove that such quantities display the properties 



E 

kik2k3 



J\^kik2k3 'pkik2k3* _ 



_|_ ^ J\^kik2k3 'j^kik2k3* 
ki k2 k3 

_|_ 'y ^ -pfclfc2fe3 J^kik2k3* 
ki k2 k3 

_ y ^ 'pkik2k3 'j^kik2k3* 
ki k2 k3 

_|_ y ^ ■^fclfc2fc3 ^kik2 

ki k2 k3 
kik2k3 

(73b) 



2^3' 



lfc2fc3 'pklk2k3* 

' b 



Similarly, for A^, Q and 5 one has 



J^kik2k3* j^ki 

kik2k3 



^2^3 



^ ' Qkik2k3* Qki 

kik2k3 
^ ^ gkik2k3* gki 

k\k2k3 



^2^3 



fe2fc3 



(74a) 



and 



_^felfe2fe3* Qk 



kik2k3 



_l_ ^ ^ J\J'k\k2k3* ^k\k2k3 
kik2 k3 

y ^ Qkik2k3* j^kik2k3 
kik2k3 

_ y ^ Qkik2k3* gkik2k3 
kik2k3 

^kik2k3* j^kik2k3 

k\k2k3 

y ^ Qkik2k3* Qkik2k3 



kik2 ks 



(74b) 



Ej^kik2k3 J^k\k2k3* _ _|_ \ ^ •pfeife2fc3 •pfclfc2 
kik2k3 kik2k3 

= + E ^, 

^1 ^2^3 



fc3'* 



fel fe2fc3 ' 
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(73a) as 



Analogous properties can be derived for terms mixing 
{M,r,n} and {Af, Q,S}. 

Let us now consider E^^, whose second-order contribu- 
tions, evaluated in Eqs. (jB28P and ()B30p . can be written 
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-11(2')/, ,N _ 



J^kik2k3 ^^fcife2fc3^* (J^kik2k3-^* j^kik^k-i 



2 i.t1 ^ - Ek,k2k3 + "iV ^ + Ek^k2k3 - i-q 

kik2k3 \ 



-.11(2") 



- Ekik2k3 + ill oj + Ekik2k3 - 



(75) 
(76) 



where the notation Ek^k2k3 = '^fci +'^fe2 +'^fe3 h^-s been introduced. Summing the two terms and using properties (j73p 
and (|74p one obtains 



Qkik2k3 ^Qkik2k3\^* ^q-)kik2k3^* jykik2k3 



kik2k3 



io - Ek^k2k3 + ill OJ + Ekj^k2k3 - iV J 



where 



Qk-ik2k3 _ r^^/cifc2^3 Jykik2k3 ^^1^2^31 

a ^ a 

q-yk\k2k3 _ ^ \\l'k-ik2k3 Qk\k2k3 ^A:iA,'2A:3l 

a ^ a 



(77) 

(78a) 
(78b) 



Notice that from Eqs. ^ and ^ follow C^i^^fcs = cfci'=2fe3 j^j-^^j 2?^ife2fe3 ^ _^^j)kik2k3 ^ ^^.j^gj. second-order 
self-energies computed in Section [B 21 can be written similarly according to 



^1 ^2^3 



^T^'\-)= E 



k\k2k3 



.22(2'+2") 



H= E 

^1^'2^'3 



CJ - 


- Ek^k2k3 +111 


" — — 

w -1- Ekik2k3 - T-V 






^^fci/c2A;3^* ^k\k2k3 


a; - 


Ekik2k3 +iv 


~^ LU + Ek^k2k3 - iV 




^'2 ^'3 ^2^^1^2^3j* 


( pkxkok^ \* nk\k2k2 
1 I'-a " ) 



oj - Ekik2k3 + iv uj + Ekik2k3 - ir] 



(79a) 
(79b) 
(79c) 



E. Energy-independent eigenvalue problem 

Defining quantities W and Z through 

(^fc - Ek,k2k3)Wt'''''- EE ^ [(C^'=^'^-^)*Z^,^ + (P^^i^-^'^.^)* v^fc] 

a 

(Cfc i?fc,fe2fc3) 2,'^'=^'=^ ^ E K""'"' + ^'''"'^ ' 

a 

Gorkov's equations f)6ip computed at second-order can be rewritten as 



^fc^a = E [(^-b ~l^Sab + Kb)Ut + Kb V^] + \cl^k2k3 y^k.k2k3 ^ (^^1/02/03)* 2,^lfc2fc3 



fcife2fc3 



^k = ^ [hl,Ut - (Tab - ^iSab + a;-,) V,'] + ^ [2?^'=^'^-^ Wt^''^''^ + (C^'=^'=■^)* Z^'^-^'^^ 

b kik2k3 



(80a) 
(80b) 

(81a) 
(81b) 



The four relations above provide a set of coupled equations for unknowns iY, V, W and Z that can be displayed in a 
matrix form 



V 

w 



+ A 




c 


/it 


-r + ;u- A* 


V 


ct 




E 










C V*\ / u 











(a- 









(82) 



r 



where is an energy-independent Hermitian matrix. The diagonalization of S is equivalent to solving Gorkov's 
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equation. Such a transformation is made possible 
by the explicit energy dependence embodied in the 
Lehmann representation, i.e. the known pole structure 
of the propagators-and consequently of second-order self- 
energies-is used to recast Gorkov's equations under the 
form of an energy-independent eigenvalue problem whose 
eigenvalues and eigenvectors yield the complete set of 
poles of Gorkov-Green's functions. The solution of such 
an eigenvalue problem has to be found self-consistently 
while satisfying Eq. 

A normalization condition for the column vectors in 
Eq. (|82p is obtained for each solution k by inserting 
second-order self-energies ([77)) and ([79]) into Eq. (|64a|) 
(or equivalently into Eq. ()64b[) ). One obtains 



E 



-E 

ki k2 ks 



\yyk1k2k3 



^kik2k3 
^k 



= 1. 



^ (83) 

The fact that ^ is Hermitian implies that eigenvalues Wfc 
are real. Moreover, similarly to the HFB problem [sg} . 
solutions come in pairs with opposite sign, i.e. for any so- 
lution {U'' ,V'' ,Wk, 2k,u!k} there exists another solution 
|-pfc*^^fc» 2*, W*, -Wfc}. This can be checked either by 
substituting w with — w in the steps that led to Eq. ([5^ 
or by re-deriving Eq. ([82)) starting from Eq. ([62]) instead 
ofEq. ([CT]). 

Let us discuss in some detail the structure of S. The 
upper-left block 



,(1) ^( T-^- 



A 



-T 



h 



A* 



(84) 



represents the "mean-field" sector. If second-order self- 
energies are zero, S = S'^-' = S^^^ and one recovers the 
HFB eigenvalue problem of fixed dimensionality (twice 
the size of the single- particle basis) for amplitudes U and 
V discussed in Sec. IIV CI 
The upper-right 



C V* 

V c* 



(85) 



prevent the achievement of convergence in an actual cal- 
culation. In practice, one limits the growth of the number 
of poles by Krylov projection techniques [13, HI], as dis- 
cussed in Paper II, while ensuring the convergence of the 
calculation. 



F. Application to = 0^ states 

The results obtained so far are general and valid for 
any choice of single-particle basis {a^} and even-number 
parity state |5'o)- If the target system, however, possesses 
specific symmetries, one can exploit them to simplify the 
set of equations. The first applications of the scheme de- 
veloped in the present paper will be dedicated to study- 
ing the ground-state of even-even semi-magic nuclei, i.e. 
states characterized by angular momentum and parity 
= 0+. Appendix [C] specifies the set of equations 
provided above to the particular case of such = 0+ 
states. 



V. QUANTITIES OF INTEREST 

A. Binding energy 

The energy sum rule first derived by Galitskii [s^ and 
formalized by Koltun j60j expresses the expectation value 
of the Hamiltonian in terms of the one-body propaga- 
tor. It is one of the appealing features of Green's func- 
tions theory, since the energy of the system, a two-body 
observable, can be computed exactly from a one-body 
quantity. 

The purpose of this subsection is to derive the analo- 
gous of the Koltun sum rule in the more general context 
of Gorkov-Green's functions. Let us first recall that the 
equation of motion of annihilation operators defined in 
their Heisenberg representation through Eq. ([23| reads 



■ daa{t) . , 



(86) 



and lower-left S^^^ blocks contain second-order couplings 
between one quasi-particlc and three-quasi-particle con- 
figurations. Such couplings further fragment the single- 
particle strength as compared to the pure HFB approxi- 
mation. As a matter of fact, following the iterative pro- 
cess leading to a self-consistent solution of Gorkov's equa- 
tions, the dimension of S^^^ grows, i.e. a larger number 
of poles is generated in Gorkov-Green's functions at each 
iteration. A propagator with an initial number of poles 
Np =2Nb generates at first iteration a second-order self- 
energy with approximately {N^f poles, which refiects 
into a matrix S^^^ of dimension -f 2(A'';,)^. After n it- 
erations the propagator and the second-order self-energy 
wiU contain respectively ©((A/'f,)^" ) and ©((A^b)^*"^^' ) 

poles, and the dimension of S'^^ will be ©((A^b)^' ^ ')■ 
This exponential growth of the number of poles seems to Also, it is useful for the following to compute the three 



and that the normal Gorkov propagator is defined at 
equal times through 

(vl/o|4(0)a,(0|*o) dujGiliio). (87) 

t=o 27ri Jct 



From the definition of the Fourier transform one can then 
derive 



dGllit) 



dt 



2711 Jc-t 
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commutators 



al [T, Ga] = -T , 

a 

a 

^4 [N,aa]=-N 

a 

Let us now write fio = (5'o|f^|5'o) as 



NN 



(89a) 
(89b) 
(89c) 



= ^0 = So 

Using Eq. (|87|) one has 

= X [r„b - M Sat] ^ f dcv Gil (^) , (90) 

while using Eqs. (|88p and ([5^ . one can also write 

Bo = -J2{'^o\al[n,aa]\^o) 

a 

daa{t) 



t=0 



(91) 



Hence one obtains the generalized Koltun sum rule 

f^O = V / diJ Gil (^) [Tab -t^6ab+UJ Sab] , (92) 

^ 47ri Jet 



where the normal Gorkov propagator G^^ appears in- 
stead of the Dyson one. 



B. Effective single-particle energies 

In Ref. [ill, an extensive discussion about effective 
single-particle energies (ESPE) in doubly-closed shell nu- 
clei was proposed. Results were based on the definition 
of ESPE going back to Baranger [6l| and the fact that 
eigenstates of the nuclear Hamiltonian are also eigen- 
states of the particle number operator. Such a defini- 
tion and its associated properties need to be revisited 
in the context of Gorkov-Green's function were particle- 
number, good symmetry, is lost, i.e. for methods 
formulated over Fock space rather than over the Hilbert 
space associated with a definite number of particles. 

As in Refs. [H, |6l|, ESPE are naturally computed as 
eigenvalues of the so-called centroid matrix, which in the 
present context is nothing but the normal part of the first 
moment of the spectral function introduced in Eq. (j55p 

^4b " ^ E l^b* E+ + Vf E~ . 



(93) 



By definition of M^^^ the computation (or extraction) 
of ESPEs requires the full spectroscopic strength, i.e. the 
complete set of separation energies and spectroscopic am- 
plitudes from both one-nucleon stripping and pickup pro- 
cesses. This is particularly critical as one moves away 
from doubly closed-shell nuclei as the low-lying strength 
becomes more and more fragmented, e.g. by pairing cor- 
relations, into both the additional and the removal chan- 
nels. This is precisely the focus of the presently developed 
theoretical method to access the complete spectroscopic 
one-nucleon addition and removal strength in open-shell 
nuclei from which ESPEs can be extracted. 

Let us now derive a sum rule for M^^-* that com- 
plements the one provided for M*^°' in Eq. (|56p and 
that provides ESPEs with a transparent physical mean- 
ing. Considering the first term on the right-hand side 
of Eq. (|93|) . substituting the definition of spectroscopic 
amplitudes and one-nucleon additional energy, one finds 

fc fc 

= E(*o|aaf^fc|«'fe)(«'fc|ail«'o) 

k 

-E(*o|f^oaa|*fe)(«'fc|aj;i*o) 

k 

-hE/^(*o|aa|*fe)(*fcl4|*o) 
fe 

^ (*o|[aa, f^]aj;i*o) + (*o|aa4l*o) , (94) 

where a completeness relation over T was removed. Sim- 
ilarly, one obtains 

Y,V^a*^bE^ = (*o|alK,f^]|*o)+M(*o|aIaa|*o),(95) 

k 

which, combined with Eqs. ([^5)) and (|94p . leads to 

m'^^^^ = {■^o]{[aaM,al}]^o) + l^Sab ■ (96) 

Using the second quantized form of T, N and , to- 
gether with symmetries of interaction matrix elements, 
one eventually obtains the key result 



^^^ab — -'-ab + 2^ Vadbc Pcd 



cd 



h°° 



(97a) 
(97b) 



where /i^ = Tab + involves the energy-independent 
(or static) part of the normal self-energy. Eventually, 
solving 



ESPE 



,ESPE „;,ESPE 



(98) 



provides ESPEs and associated single-particle wave func- 
tions. Defined in this way, ESPEs are manifestly in- 
dependent of the single-particle basis used to compute 
the centroid matrix (|93p . They possess the meaning of 
an average of observable one-nucleon separation energies 
weighted by the probability to reach the corresponding 
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many-body state of the iV ± 1-body system by adding or 
removing a nucleon in the single-particle state ip^^^^. As 
such, it is however essential to understand that ESPEs 
are by essence non observable and display an intrinsic 
resolution scale dependence, just as spectroscopic factors 

do nil. 

In spite of the breaking of particle-number symmetry 
leading to the coupling of additional and removal spectro- 
scopic amplitudes via anomalous self-energies in Gorkov's 
equations, Eq. (|97bp demonstrates that the centroid ma- 
trix is equal to the normal static field h°°, exactl y as for 
theories that explicitly conserve particle number [53 ■ In 
other words, the centroid sum rule does not only screen 
out the energy-dependent part of the normal self-energy 
but also screens out the entire anomalous self-energy. 
This is an a priori non trivial, though straightforward 
to obtain, result. Of course, the explicit tackling of pair- 
ing correlations through anomalous propagators and self- 
energies does impact the results indirectly via their feed- 
back onto the normal one-body density matrix p entering 

Eventually, one can demonstrate, just as for particle- 
number conserving theorie s 1551 t hat h°° only involves 
the so-called monopole part |62l . l63j of the NN interaction 
whenever |^o) is a J'^ = 0"*" state. 



C. Natural basis 

The natural basis is the one that diagonalizes the one- 
body density matrix, i.e. 

pV'r=«r^r' (99) 

where basis states i/'p^' are called "natural orbitals" and 
where diagonal elements rip^' denote "natural occupation 
numbers" . The A'' most occupied natural orbitals define 
the set that better approximates the true (correlated) 
density matrix p in terms of a Slater determinant wave 
function. Thus, the natural basis is most convenient to 
expand approximations of obscrvablcs other than ener- 
gies (e.g., radii, density distributions, etc). We stress 
that in general i/jp*^* correspond to superpositions of or- 
bits V'p with ESPEs both above and below the Fermi 
surface, chosen to optimize the density profile of the sys- 
tem, and therefore are a poorer approximation to energy 
levels. Conversely, single-particle states ip^^^'^, in Eq. 
((98)) , can be directly associated to orbits of the effective 
shell structure. 



D. One-body observables and radii 

The expectation value of a general one-body opera- 
tor, O, is obtained from the normal density matrix (|54a|) 
or, equivalcntly, by integration over the normal Gorkov 



propagator 

ah -^Ct 

(100) 

where Oah — {a\0\b) are matrix elements of the one-body 
operator. 

For matter radii, however, one needs to sum over par- 
ticle positions in the intrinsic frame. The operator for 
the root mean square point radius rrms is thus 

^rms — ^ ^ R-) 

i 

= l('l_lWr2 - J-Vr^.r,, (101) 

where denote coordinates of nucleon i in the labora- 
tory frame and R = {J^i ^i)/^ is the centre-of-mass co- 
ordinate. Operator (jlOip contains a two-body correction 
term and depends on the number of particles. As for the 
centre-of-mass corrected Hamiltonian, this is the form 
suitable for applications in Fock space. At second or- 
der in the Gorkov self-energy, as considered in this work, 
there is no resummation of diagrams corresponding to 
correlated two-body propagators. Only the free propa- 
gation of dressed pp or hh is accounted for in the second 
order diagrams. Correspondingly, one can approximate 
the two-particle density matrix with the antisymmetrizcd 
product of correlated one-body density matrices in order 
to evaluate the two-body part in (|10ip . The rrms radius 
is thus calculated to first order in A^^^ as 

where r^j, = (ajr^jb) and where {ab\ri-r2\cd) are antisym- 
metrizcd two-body matrix elements while N denotes the 
total number of nucleons (protons plus neutrons). Eqs. 
(jlOip and (|102p are valid for matter radii. Charge point 
radii r^h can be obtained from 

^ p 

where Z is the number operator for protons and where p 
runs only over protons while R is the same as in (jlOip . 
Isotope shifts are calculated from differences of squared 
charge radii, 

Sirl.r^''' ^ {rl,r' - {rl,r , (104) 

where A' is the number of nucleons of the system un- 
der consideration whereas A^ characterizes a reference 
nucleus. 



16 



E. Pairing gaps 

Experimentally, a suitable way of extracting the pair- 
ing gap goes through, e.g., the three-point mass formula 

A(3)(iV) = - 2E^ + , (105) 

where N is the total number of nucleons. This is moti- 
vated by the relation between the odd-even staggering of 
nuclear binding energies and the lack of binding of the 
unpaired odd nucleon, as first pointed out in Ref. [6^ . 

To compute A*^'^^(7V) theoretically, one needs to per- 
form consistent calculations of odd nuclei. In the present 
context, this would require to perform Gorkov calcula- 
tions for a state |^'o) having an odd number-parity quan- 
tum number, i.e. a state such that the sum runs over 
odd N in Eq. p6|) . This is however beyond the scope of 
the present work. The next best approximation would 
consist in keeping an even number-parity state while ac- 
counting for the blocking of a quasi-particle within the 
filling approximation |65|- Such an approximation re- 
mains however to be formulated within the general frame 
of Gorkov-Green's function formalism. 

In such a situation, the next best estimate to the 
ground-state energy of the odd system is obtained 
through [H, [63 

(106) 

where * is the energy of the odd nucleus computed as 
it were an even one, i.e. as a fully paired vacuum with 
an odd number of particles on average, while cup denotes 
the lowest pole energy obtained from that Gorkov calcu- 
lation. Obviously, for even N one simply has Eq = Eq * . 
With such an appropriate decomposition of the energy, 
the three-point mass formula reads 

A(3)(iV)«L_^ ^ + A^(iV). (107) 

The first contribution relates to the second derivative of 
the smooth part of the energy E^*, i.e. the energy 
curve on which both even and odd nuclei would lie in 
the absence of odd-even mass staggering. Such a sec- 
ond derivative of Eq * evolves very smoothly with N. 
However, the corresponding contribution to A^'^^ (N) os- 
cillates strongly around zero due to the factor (—1)^ ap- 
pearing in Eq. (jl07p . accounting for the odd-even oscil- 
lation of experimental A^'^^ (N) and having nothing to 
do with the pairing gap itself [i^ H^l • The second con- 
tribution to A^'^)(A^) relates specifically to the unpaired 
character of the odd nucleon, and thus extracts the ac- 
tual pairing gap at the Fermi energy in open-shell nu- 
clei [ei 113 

A /,rN f for odd 

[(c^p +ujp^ )/2 for A even 



VI. $-DERIVABILITY 

Gorkov-Green's functions constitute a versatile and 
powerful technique that can as well be applied to time- 
dependent, non-equilibrium systems. When truncating 
the self-energy expansion, and hence approximating the 
solution of the many-problem, one has to pay attention 
to the possible violation of basic conservation laws in- 
volving e.g. particle number, total energy, total momen- 
tum, total angular momentum. For an arbitrary set of 
self-energy diagrams nothing assures that quantities are 
conserved with time even though the corresponding op- 
erator commutes with H. 

A way to construct a class of conserving approxi- 
mations was devised in the early 1960's by Baym and 
Kadanoff [6^, H^. Baym and Kadanoff demonstrated 
that if the self-energy is derived from a certain functional 
$ of the one-body Green's function (and the two-body 
potential) , previously introduced by Luttinger and Ward 
[70| , the resulting scheme automatically satisfies all basic 
conservation laws. Moreover, the resulting approxima- 
tion preserves thermodynamic consistency requirements, 
including the Hugenhotlz-van Hove [7l| and Luttinger 
[t^ identities. Fulfilling such consistency requirements 
avoids ambiguities in the calculation of thermodynamic 
observablcs, i.e. different ways of computing the same 
quantity yield the same result. 

The concept of ^-derivable approximations, intro- 
duced by Baym and Kadanoff for normal Green's func- 
tions, was generalized to Gorkov's formalism by De Do- 
minicis and Martin [t^ (see also [zl, [t^ ) . In this case 
such a class of approximations relies on the existence of 
a closed functional $ of the four Gorkov-Green's func- 
tions G and the two-body interaction 1/, from which 
self-energy contributions are obtained via a functional 
derivative. At zero temperature, such a functional is 
closely related to the correlation energy AEq, i.e. the 
total e nerg y measured with respect to the unperturbed 
energy [7y| 

A£;o = $[G,F]+Tr{§^-l}-Tr{ln|^}, 

(109) 

where traces have to be performed over Gorkov and one- 
particle Hilbert spaces 

Tr.Tr«,iTr,/^. (110) 

There exist two possible strategies to build a $- 
derivable scheme. Starting from a carefully chosen set 
of self-energy contributions, Baym-Kadanoff's functional 
can be formally defined through [76j 

<i>[G,F]^-- 5^-Tr{s(")[G,F]G} 

OO 

= ^$(")[G,X/], (111) 
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where S''"^[G,1^] denotes skeleton self-energy terms of 
order n. Compared to the standard definition, an addi- 
tional 1/2 factor appears in relation to the trace over the 
two-dimensional Gorkov space in Eq. (jllOp . Notice that 
self-consistency, i.e. the use of dressed propagators in the 
functional, is a necessary condition for $-derivability. 

Alternatively, one can use diagrammatic techniques to 
construct analogously to the self-energy expansion. 
At order n in 1/, is given by two-fermion-line irre- 
ducible connected closed skeleton diagrams. With the 
obvious change from an open to a closed topology, all di- 
agrammatic rules outlined in Appendix IB II hold for the 
construction of the functional. From $ one can obtain 
the four self-energies by differentiating with respect to 
Gorkov propagators 



(5$[G,V] 



(112a) 









Sl^H = -2 



(112b) 



(112c) 



(112d) 



FIG. 5. Diagrams contributing to $^^K 



that generate first- and second-order sclf-encrgics of Figs. 
[U [5] and [3] are depicted respectively in Figs. U and [S] 
Writing explicit expressions for such diagrams and ap- 
plying functional derivative (|112|) . one obtains all self- 
energies calculated in Appendix IB 21 As an illustration, 
the first order case is treated in full details in Appendix 
ID2] 



as demonstrated in Appendix ID II Any subset of $ dia- 
grams employed to derive the self-energy via Eq. (|112p 
will generate a conserving approximation. If at a given 
order all terms are taken into account, the resulting self- 
energy will contain all possible contributions at that or- 
der. Eventually, one should have internal consistency, 
i.e. if the functional is constructed from a certain set 
of self-energy contributions via Eq. (|llll) . all and only 
these self-energy contributions must be generated from 
that functional when employing Eq. (|112p . 






FIG. 4. Diagrams contributing to <I>'^\ 

Most of the commonly used (fully) self-consistent ap- 
proximations in Green's function theory are ^-derivable. 
It is the case of the second-order approximation used in 
the present work, as well as of Hartree-Fock, RPA or 
T-matrix approximations. Diagrams in the ^-functional 



VII. CONCLUSIONS 

The aim of the present paper is to extend the reach 
of ab-initio nuclear structure calculations to truly open 
shell nuclei. This is done by implementing self-consistent 
Green's Function method within the general Gorkov's 
scheme. Such a method retains the simplicity of single- 
reference approaches, i.e. methods relying on a many- 
body expansion around a single vacuum. As for open- 
shell systems, this can only be done at the price of break- 
ing C/(l) symmetry associated with particle- number con- 
servation in order to tackle Cooper pair instabilities and 
explicitly account for pairing correlations. 

The present work provides a detailed account of 
Gorkov's formalism and of its specification to second- 
order in the expansion of normal and anomalous self- 
energies. At the present stage, this is done in terms 
of two-nuclcon interactions only. First numerical ap- 
plications of such a scheme will be reported soon in a 
forthcoming publication. This constitutes the first ab- 
initio application of self-consistent Gorkov-Green's func- 
tion method in finite nuclei using realistic two-nucleon 
interactions. The extension of the proposed ab-initio 
method to more advanced truncation schemes and to in- 
clude threc-nuclcon forces is the aim of future works. 



18 



ACKNOWLEDGEMENTS 



with 



f{z)dz 



(A2) 



The authors would hke to thank S. Baroni and P. 1 
Bozck for useful discussions. This work was supported ^" ~ 27ri 

by the United Kingdom Science and Technology Facilities 
Council (STFC) through travel grant No. ST/1003363. 

V. S. acknowledges support from Espace de Structure and where C is a contour contammg c on which / is an 



Nucleairc Thcorique (ESNT). 



Appendix A: Normalization condition 

In order to work out the normalization of the spec- 
troscopic amplitudes let us consider the expansion of 
Gorkov's equation ((34)) around the pole +ijJk- We re- 
call that a complex function f{z) can be expanded in a 
Laurent series around a point c in the complex plane as 



-t-oo 



f{z) = X! ^n{z- CY 



(Al) 



alytic. Here one is interested in the case in which f{z) 
has a simple pole at z ^ c, which means that the inte- 
grand in Eq. (|A2[) has a pole of order n + 2 at z = c and 
the integral can be performed by means of the residue 
theorem. In particular, the ti = coefficient reads 



ao 



1 



f(z)dz d 



(z-c) 



2 fjz) 

{z^c) 



(A3) 



Performing a Laurent expansion of Gorkov's equation 
P4p around tu = -l-Wfc and extracting the coefficients of 
order zero on both sides, one obtains 



lim — 



cd 



+(^ - ujk) - nu)-', Kbi^)] } , (A4) 



where the singular part in the dressed propagator has been isolated as follows 



■vk' Y*:'t 

uj -ujk + irj 'j-f LO - LOk' -[- i 



17) 



Applying the derivative to all terms and dropping the ones that give zero in the limit a; — > Wfe, one has 

cd 



(A5) 



hm {B^,(..) = - ^u)ab + E f - - ^vTa! ^^'^(^) ^^^'^ 



duj 



■ X^Xf + _ ^^-^-1 s,,(c.) } , (A6) 



such that, using Eq. (|6T|) . 



lim { E - ^u);a X^X^-t + = (c. - Qu);^ + E f " 

^— ^^fc — — L 



1 d'Scdji^) 



(A7) 



Multiplying both sides of Eq. (|A7[) by [uj — ^u)fe and where the terms involving B*^ cancel out after using the 



summing over e yields 



(AS) 



conjugate Gorkov's equation 
Gabico) = Gi°)(c.) + E Gacico) t^^^) G^"^ (c.) . (A9) 



cd 



Multiplying by X^^ from the left, summing over / and 
renaming (/, d) to (a, b) one finally obtains the normal- 
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ization condition 



second to the first index 



a ab 



Xk 
b ' 



(AlO) 



where only the proper self-energy appears as a result of 
the energy independence of the auxiliary potential. Sim- 
ilarly one can derive a condition for Gorkov's amplitude 
Y 



, (B2a) 



E 



Y„N =1 



ab 



(All) 



, (B2b) 



Appendix B: Diagrammatic 

1. Diagrammatic rules 

A convenient way to express the expansion of the 
single-particle propagator is via diagrammatic tech- 
niques. By giving the interaction and the single-particle 
propagator a graphical representation and by establish- 
ing a set of rules one can generate diagrams that are in 
one-to-one correspondence with the terms appearing in 
the expansion. As it provides an immediate insight to 
physical processes associated with the various contribu- 
tions, the diagrammatic expansion is of great help when 
choosing a suitable approximation. It is relevant to dis- 
cuss diagrammatic rules in some details here given that 
there exist differences compared to rules applicable to 
the diagrammatic expansion involving normal contrac- 
tions only. 

In the present work antisymmetrizcd interaction ma- 
trix elements are represented by a dashed line labeled by 
four single-particle indices 



a b 

Vabcd = • • . (Bl) 

C d 



a 

Gll^"\co)^ I , (B2c) 

b 

a 

gT\uj)^ I . (B2d) 

b 



One should notice that, as opposed to traditional graph- 
ical representations of Dyson's propagator, Gorkov's 
propagators carry two arrows specifying whether a given 
propagator results from the contraction of two creation 
operators, of two annihilation operators, or of one cre- 
ation (annihilation) and one annihilation (creation) op- 
erator. 

With building blocks (IBip and (|B2p one can construct, 
order by order, the (diagrammatic) perturbative expan- 
sion for each of the four Gorkov propagators ((22)) . To 
obtain all terms of the expansion at a certain order m 
and for one of the four Gorkov propagators, the follow- 
ing rules are employed: 



Single-particle unperturbed propagators, i.e. Green's 
functions associated with the unperturbed Hamiltonian 
ill/ introduced in Eq. (PT|) . are depicted as solid lines 
labelled by two indices and one energy flowing from the 



1. Draw all topologically distinct connected direct di- 
agrams with m horizontal interaction lines (with 4 
single-particle indices) and 2m + 1 directed prop- 
agation lines (with 2 single-particle indices each, 
connecting the 4to indices of the interaction and 
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the 2 external ones) 



( 

\ t 



(B3a) 



(B3b) 




(B6) 



Notice that exactly two incoming and two outgoing 
lines must be attached to a given interaction vertex, 
i.e. diagram (|B3a[) is allowed while diagram (jB3b[) 
is not. 

Topologically distinct diagrams cannot be trans- 
formed into each other by any translation (in 
the two-dimensional plane) of any of the vertices 
without disconnecting or reconnecting propagation 
lines. 




(B4) 




(B5) 



For example second-order diagrams (jB4p and (jBSP 
are topologically equivalent, while diagram (jB6[) is 
not. 



Connected diagrams are diagrams in which it is 
possible to go from each interaction line to any 
other by moving along propagation lines. 




(B7) 



1 



(B8) 



For example first-order diagram (|B7|) is connected 
while first-order diagram (|B8|) is disconnected. 



For a given diagram, exchange diagrams are derived 
by exchanging the end points of two propagation 
lines coming in or out of one or more interaction 
vertices. Since we are using anti symmetrized ma- 
trix elements in (|Bip , it follows that for each set of 
diagrams obtainable from one another by means of 
such exchanges, one must only retain one represen- 
tative diagram, arbitrarily chosen and denoted as 
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direct, and discard all the other ones. 



diagram. 




(B9) 



For example if one considers diagram (|B6p as direct 
(the choice of the present work) it follows that one 
must discard diagram (|B9[) . 

In cases where it is unclear whether diagrams are 
topologically distinct, one can always resort to a 
direct application of Wick's theorem. 

2. Assign an energy to all propagation lines such that 
the energy in each interaction is conserved (the en- 
ergy entering a vertex must be equal to the energy 
exiting). As a result, a m-order diagram will have 
m independent internal energies and the incoming 
external energy will be equal to the outgoing exter- 
nal one. For each independent energy, multiply by 
a factor l/27r. 

3. Write down a V (with corresponding s.p. indices) 
for each interaction line and a G^^^^ (with corre- 
sponding s.p. indices and energy) for each prop- 
agation linqj according to representation ()B2p . If 
the energy tu flowing along the propagator has the 
opposite direction than in definition (|B2p . the as- 
sociated term is G^^^^{-~ijj). 

4. Write an overall factor i™. 

5. Write a factor 1/2 for each pair of equivalent prop- 
agation lines, i.e. pairs of lines starting at the same 
interaction vertex and ending at the same interac- 
tion vertex and corresponding to the same type of 
Gorkov propagator. This factor is due to the an- 
tisymmetrization of the potential, i.e. to the fact 
that exchanging the incoming lines of two interac- 
tions connected by equivalent lines yields the same 



1 

I 

( I • 

;; 

( I ' 




(BIO) 



For example diagram (jB6[) has a pair of equivalent 
lines, i.e. those labeled by (c, J) and {f,h), while 
diagram (jBlOp has none. 

6. Write a factor 1/2 for each anomalous propagator 
starting and ending at the same interaction ver- 
tex. This factor appears for the reason discussed in 
point 5 and applies, e.g., to diagram (|B7p . 

7. Write a factor (— l)^c+^» where Nc is the number 
of closed fermionic loops and Na is the number of 
anomalous contractions. 

8. Interpret equal-time propagators as 

\jmGllit,t')=GiliO,^rj), (Blla) 

]imGll{t,t')=Gll{0,-rj), (Bllb) 

= G2i(0,-r?), (Bile) 

= Gll{0,+r^), (Blld) 



limGllit,t') 
\jrnGlUt,t') 



which implies that integrations over uj are per- 
formed in the complex energy plane, either by clos- 
ing the contour in the upper {C t) or in the lower 
(C l) half plane as 



1 dcuGiliiu)- 


Jet 


dujGiliu;), 


(B12a) 


1 dcuGlliiu)- 


Jet 


dioGlliio), 


(B12b) 


1 dcoGlliu)- 


Jet 


duGlliu), 


(B12c) 


1 duGlliu)- 


let 


dcoGllico). 


(B12d) 



Any normal propagtion line can be interpreted either as G^^ 
or G^^ . The choice of the present work consists in identifying 
normal lines with G^^ in the expansions of G^^ , G^^ and G^^ 
while using G^^ in the expansion of G^^ . 



When equal-time propagators appear the ordering 
of the annihilation and creation operators must be 
as in the starting Hamiltonian. Hence limits (jBlla[) 
and (jBlldl) must be taken in opposite ways. Since 
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the operators in G^^ and G^^ anticommute the re- 
maining two hniits can be arbitrarily interpreted, 
as long as they arc taken consistently. 



9. Sum over all internal single-particle indices and in- 
tegrate over all internal energies. External indices 
and energy refer to the Gorkov propagator being 
expanded. 



Once the expansions of the four one-body Gorkov prop- 
agators are written down, one can derive the correspond- 
ing expansions for the self-energies by simply stripping 
off external propagation lines, e.g. to the term (jB13a[) 
corresponds self-energy contribution (|B13b[) . 



t w a 



i 



tw b 



(B13a) 



(B13b) 



d J W/ 



--0 



loj' (B14b) 



Once this distinction is made, one can demonstrate 
that the complete propagator expansion is generated by 
keeping irreducible skeleton self-energy diagrams only 
and by replacing accordingly in such diagrams all un- 
perturbed propagators by dressed ones. Dressed propa- 
gators are Green's functions that are solution of Gorkov's 
equations: their appearance in the self-energy expansion 
generates the self-consistency characterizing the method. 

It thus follows that only irreducible skeleton self-energy 
diagrams with dressed propagators have to be computed. 
Such dressed propagators are depicted as solid double 
lines and are labeled by two indices as well as by an 
energy, just as for unperturbed ones, i.e. 



G 



(B15a) 



All self-energy contributions can be divided into two 
types: one-line reducible and irreducible self-energies. Ir- 
reducible self-energies are constituted by diagrams that 
can not be separated into two parts by cutting one propa- 
gation line. For example diagram (jB4p is reducible while 
diagram (|B6[) is irreducible. Irreducible contributions 
can be further divided into skeleton and composed di- 
agrams. Skeleton (composed) self-energies are obtained 
by keeping, at a given order order, only those terms that 
cannot (can) be generated by successive insertions of ir- 
reducible self-energy terms of lower order. At first order, 
all diagrams are irreducible by definition. An example 
at second order is given by the two diagrams (|B14[) : the 
first term (|B14a|) is a skeleton diagram while the second 
self-energy contribution (jB14b[) can be generated by two 
successive insertions of the first-order term ()B13bp . 



, (B15b) 



, (B15c) 



(B15d) 




(B14a) 



Diagrammatic rules to compute irreducible self-energies 
are the same as for reducible ones, with the only dif- 
ference that dressed propagators (|B15[) have to be used 
instead of unperturbed ones. 
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2. Self-energies 



which reads 



The present section addresses the derivation of first- 
and second-order self-energy diagrams. 



a. First order 



The first normal contribution corresponds to the stan- 
dard Hartree-Fock self-energy. It is depicted as 



.11(1) 



a c 
• 

h d 




(B16) 



and reads 



.11(1) 



I ^J2^acMGYA^'), (B17) 



where the energy integral is to be performed in the up- 
per half of the complex energy plane, according to the 
convention introduced in Rule 8. Inserting the Lehmann 
form pSap of the propagator one obtains 



-22(1)^, A _ 



''ab 



(uj) = -i 



o / J ' oaac I I 



cd^k 

^_5.ii(i) 

oa 



(B20) 



The anomalous contributions to the self-energy at first 
order are 



.12(1) 
-'ab 



(^) 




^ LU' 



, (B21) 



duj' X - - 

"-'T cd,k 



■bd 



- « / Vacbd ■ , , 

let 271" ^ W' + Wfc - IT] 

cd.k 



(B18) 



.21(1) 
^ab 




, (B22) 



where the residue theorem has been used, i.e. the first 
term, with +if] in the denominator, contains no pole in 
the upper plane and thus cancels out. As in the standard 
case the Hartree-Fock self-energy is energy independent. 

Similarly one computes the other normal self-energy 
term 



.22 (1) 



a c 
• 

b d 




(B19) 



and are written respectively as 



2 Jet 27r ^ 

o / o / J *^ abed I I 

2 Jet 27r ^ w' - ujk + IV 



(B23) 



cd.k 
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and 



been used. 



,21(1) 
-•ab 



Veda!, Gldi^') 

'-'T cd 



duj' 



r, / , ^cdab ■ , . 



duj' 



uj' -\- iUk — if] 



cd,k 



cd,k 

^12(l)i* 



h. Second order 

Let US now proceed to the computation of the second- 
order contributions. The first term is the standard 
second-order self-energy 



= (B24) 

where the same integration technique as in (jB18| has reading 
I 



d 
b 




iuj'" (B25) 



T'\^) = \j E VaecfVd,bHGl\{u')Gf^{u.'')Gi\{u''')5{u,- 

cdefgh 

= W E ^-^/ ^'^^^^ Gl\{u:') Gf^iuj") Gl\{J + u" - c.) . (B26) 



Notice that the minus sign coming from rule 4 is cancelled by a minus sign coming from the presence of a closed loop 
(rule 7). The integrations over the two energy variables are performed in this case using two successive applications 
of the formula 



m = — 



2^ \E' 
F1B2 



fi + iri 



E' - hi 
F2B1 



irj 



Fo 



Bo 



E'-E-f2 



E' -E-bo 



IT] 



(/i-62) + w; E + {f2 - bi) ~ IT] 



(B27) 



The above integral, defined on the real axis, is computed by extending the integration to a large semicircle in the 
upper or lower complex half plane of E' (this can be done since the integrand behaves as for \E'\ — )- 00 and 

this branch do not contribute to the integral) and then by using the residue theorem. Of the four terms, two have 
poles in the same half plane and yield zero as the contour can be closed in the other half. Applying this formula to 
the integral (jB26l) one obtains 



^11(2'). . 1 fdu'dc" / U!:^U',^* -P^^-P^ 

caej gh,kik2k^ K. . 





IT] 



V Vaecf Vdgbh I / / \ \ , + , / / ' \ . . (B28) 



cdef gh,kik2 fca 
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With the same technique one can evaluate all other terms contributing to the second-order self-energy. One has 




(B29) 



reading 



E V.ecfV,,,j,G\^{J)Gf,{J')Gl\{.' +J' -u) (B30) 

cdef gh 



The two diagrams of the other normal self-energy are respectively 



yielding 




(B31) 
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cdef gh 



E 



cdef gh,kik2k3 
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(B32) 



and 



reading 



^22 (2" 




(B33) 



dcu' did'' 



E ^'^^ ^S/^.^-^. ^c^(^') Gfg{uj") Gl){u' + u" ^) 



cdefgh 

VcfhE VhhrlT, \ -, ^ : h 
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(B34) 
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The first of the anomalous self-energy is 




for what concerns the first contribution, which reads 



^T'\^) = -J E Kec/W-G^^(-')G)lK)Gji(a.'+..''-c.) (B36) 



cdefgh 



^aecf yhbgd 

cdefgh,kik2k3 



oj - {ujki + + ^ks ) {^k^ + ^ki + Wfe J - 177 I 



and 



sf,('")H= tc^'I t^"I Die.'" , (B37) 



for the second contribution yielding 



^T"\-)-\j E %ecfVu,,G]^A^')Gf^{u")Gl\{u:' + -u) (B38) 



cdefgh 



^fci ^fcl * ^k2 y fc2 * * y fe3 -pfel * * Z:/fc2 V*=3 ^fe3 * 



cdefgh, kik2k3 \ 



Finally 



S^;('')M= t^'D t^"ir li^'" , (B39) 



reads as 



^T'\-)--j E %f-aeVgdh»GfA^')G\]{u")Gl){u' +U:" (B40) 



E ^S/- ^a^il^b ^ '^'^ + •'d "e "g -7 
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while 



^21 (2" 




(B41) 



is expressed as 

^21(2") 
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cdefgh 
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cdefgh,kik2k^ 
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Appendix C: = 0^ states 

The present section specifies the complete set of equa- 
tions to = 0+ states. 



1. Time-reversal invariant systems 

Let us define the time-reversal operator as 

j=e^^Syj^^ (CI) 

where K is an operator that associates to a wavcfunction 
its complex conjugate and Sy is the y-axis spin-projection 
operator of the N-body system. The time-reversal oper- 
ator is antiunitary (unitary and antilinear) and displays 
the following properties 

T^T = 1 , (C2a) 

= (-1)^ , (C2b) 

{a\T)\b) = [{a\{T\b)Y , (C2c) 
(a| (rt \h) = (6| (T \a) . (C2d) 

One can also introduce the time reversal operator acting 
in Fock space. This is done by specifying its transfor- 
mation rules on standard creation and annihilation op- 
erators defined on the tensorial product of spatial (rep- 
resented by the vector r), spin (represented by the spin 
projection nic) and isospin (represented by the isospin 
projection q) one-body Hilbcrt spaces 



T'l'at r =2m^at^ , 


(C3a) 




(C3b) 


1^ drrricjq 'T —2 Tfl(j drrfi^jq i 


(C3c) 


1~ dr^n^q 7"^ —-2 fflfj (Xrrh^q i 


(C3d) 


as well as on the particle vacuum 




r|o)=|o), 


(C3e) 



I 

where the notation m„ = —m^ has been used. Given 
Eq. (|C3p it is easy to prove that for any basis {aj^} dosed 
under time-reversal, defining 



(C4) 



provides a partner basis of the type ([1]) . Accordingly we 
define the action of T on one- and two-particle Dirac ket 
and bra as 

T\a)= I a) (a|T'l' = 77a (a| , (C5a) 

T\ab) ='qa'nh\ab) (a6| T''' = 77^ 77^ (a5| . (C5b) 

It follows that for kinetic energy, which fulfills TT = 
T, 

T,,^ia\r)TiT\b) 
= [ia\{r^Tr\b)Y 
= [{a\T\b)r 

= , (C6) 

and similarly for time-reversal invariant interactions, i.e. 

Vvr^v, 

%l,i^{ab\T^)V{T{\cd)~\dc)} 
= [{ab\[T^VT{\cd)^\dc)}Y 
= [{ab\V{\cd)-\dc)}r 
= KTbcd • (C7) 
Considering a time-reversal invariant system, i.e. a refer- 
ence state \^()) satisfying T|^'o) = l^'o) and (^'o|T'^) = 
(^'ol, one can prove, using property (|C2[) . that in this par- 
ticular case the anomalous density matrix (j54cp is Her- 
mitian 

Pab = (^'olofcaal^'o) 

= (vl/oirt) (TabTt aa|*o) 

= -[{^o\a, (rt a, rtr 1*0)]* 



Plb- 



(C8) 
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2. Single-particle basis 

In the remaining of the present section the many-body 
system under study is assumed to be in a jn = 0+ 
state, where the parity is 11 = (—1)^. A possible choice 
for labehng single-particle basis states in this context 
is a = {n,£,j,m,q), where n represents the principal 
quantum number, £ is the orbital angular momentum, j 
is the total angular momentum, m is the projection of 
the total angular momentum along the z axis and q is 
the isospin projection. The spin a is omitted from the 
single-particle label because trivially ct = 1/2 for all nu- 
cleons. A choice that will appear to be more convenient 
below consists of labeling single-particle states according 
to a = {n,Tr, j,m,q), where the parity tt = ( — 1)^ sub- 
stitutes the orbital angular momentum. Since for a given 
j, £ — j ± ^ are the only possible values, there exists a 
one-to-one correspondence between £ and tt. 

Different phase conventions exist to define single- 
particle states in such a context. In the present work, 
spinors are written as 



(f |a) = ipa{rq) 
_ Unijjrq) 
r 

_ Unijirq) 



(C9) 



where C denotes a Clebsch-Gordan coefficient according 
to 



C: 



JM 

jimij2m2 



{jimij2m2\JM) 



(CIO) 
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The rijimif) are spherical spinors that recouple the an- 
gular part of the wavefunction to spin-1/2 spinors. They 
fulfill 

.P njirnif) - j{j + 1) , (Clla) 

L^nje^{r)^h^ei£ + i)n,e^if), (ciib) 

S'^^jem{f) = a{a + 1) D,jim{f) 

(Cllc) 

Jz ^jemif) = hmfljirn{f) ■ (Clld) 
Spherical spinors are orthonormal, i.e. 

27r 7r 

J d(p j M Sin(?9) ^f) ^.ytmi{'&, ^) Sjj' fe' 5ram' , 



and fulfill 



j 



J2 n],Jf)<J^},,„r{f)=0. 



(C12b) 



As shown by Eq. (jCllcp . the total spin cr is a good 
quantum number while its projection is not, since single- 
particle states mix = ±i. If the time-reversal oper- 
ator (jCip is applied to a state a belonging to this single- 
particle basis one can easily prove, using 



2m„ 



y; 



•(r) 



(C13a) 
(C13b) 



as well as standard properties of Clcbsh-Gordan coeffi- 
cients, that 

(rq) , (C14) 



where rjijm = (— 1)^"-'"'". Equation (|C14p demonstrates 
that time-reversal operation connects basis state a, up to 
a phase, to state d = (n, £, j, —to, q) and thus constitutes 
an anti-unitary transformation that can be employed to 
define the partner basis {d^} used to introduce Gorkov 
Green's functions. The creation and annihilation oper- 
ators introduced in Eq. ([1]) in this case take the form 



(C15a) 
(C15b) 



Equivalently, one obtains d = (n, tt, j, —to, q) if parity is 
chosen to label single-particle basis (which will be our 
choice in practice), such that 



'^mrjmq — VTrjm O-mrj-mq i 



(C16a) 
(C16b) 



with riT^jm = TT (— 1)-'+™. By including ry^jm in the defi- 
nition of the building blocks of the theory (density ma- 
trices, propagators, etc ...) one can exploit symmetry 
properties associated with time-reversal invariance that 
lead to a simplification of the formalism. 



3. Block-diagonal structure of propagators 

It is easy to prove that the operator a\^-^^^ {o^mrjmq) 
is the m}^ component of an irreducible tensor of rank j 
and that the corresponding annihilation operator trans- 
forms contragrediently, implying that (— l)™aj_„i is also 
the m}^ component of an irreducible tensor of rank j. 
Starting from such a property, one can demonstrate that, 
in addition to being diagonal in isospin space as only 
proton-proton and neutron-neutron pairing is considered 
here, Gorkov propagators possess a block-diagonal struc- 
ture relative to quantum numbers tt [£) and j and are 
independent on to, i.e. 



riajaiamaqaribjbibmbqb ^ ' 

_ (^9192 t \ 

UamaaribmbP^ ' 



(C17) 
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where the notation a = {ja^T^a^qa), i-e. a = (na,a,ma), 
has been introduced. Similarly, one writes unperturbed 
propagators as 



G 



9192 (0) 
ah 



,<5„^G^f,?H(o)(^). (C18) 



The Lchmann representation also reflects the block- 
diagonal form of the propagators. In particular, there ex- 
ist selection rules associated with label k = {rik, k! , mj,), 
with k' = {jk,T^k,Qk), characterizing many-body states 
l^'fc) introduced in Eq. (p6)) . Considering the definition 
of the spectroscopic amplitudes ((39|) and (|40l) and apply- 
ing Wigner-Eckart theorem, one finds 



-^00 



(-1)™° C™ 



(*o||an„a||^'n. 



V2ja + 1 



(C19a) 



where k = {ik,T^k,Qk — Qo), with Qo being the isospin 
projection of |^'o)- It is assumed here that, analogously 
to angular momentum, |^'o) has a good isospin projec- 
tion and that Qk is determined by the isospin projection 
of the creation/annihilation operator acting on |^'o), i.e. 
^Qk-Qa-qa- Froni parity conservation follows S-,:^-^^. As 
the spin of j^'fc) can be different from 1/2, the fact that 
TTfe = ttq does not imply Ik = £a, hence the single-particle 
basis is labeled by a = {na,TTa,ja,'m'a,qa) from now on. 
Similarly one has 
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^ r. 



(C19b) 



- <^JkJ. Orakra^ ^2j^ 

— 6 6 V"*^ 



:*o||aLaP«.«) 



6 6 -n- V"'' 



(C19c) 
(C19d) 



Inserting Eqs. (|C19|) into Eqs. ((SH]), the set of Gorkov 
Green's functions can be written according to Eq. (jC17[) 
as 
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W -f Wfe — 





, (C20a) 
, (C20b) 
, (C20c) 
, (C20d) 



where only one sum over principal quantum number Uk 
remains. 



4. Matrix elements of the nuclear potential 

Let us consider two-body interaction antisymmetrized matrix elements Vabcd introduced in Eq. (j4|), which depend 
on angular momenta ja-, jb, jc, jd of the two incoming and two outgoing nucleons as well as on their third components 
nia, rub, nic, md- Writing all indices explicitly they read 

= (1: {naama); 2: {ub ■mb)\V^^ \1: {uc^mc)] 2: {nd6md)) 
- (1: {nactma); 2: {nb P ■mb)\V^^ \l: {na 6 md); 2: {nc^nic)) 
= {l:{naama);2:{nbPmb)\V^^\l:{n,-im,);2:{nd6md)) . (C21) 

One can go from such a representation, referred to as the m-scheme, to the jj-coupled scheme or J-scheme, in 
which incoming and outgoing two-nucleon states are labeled by total angular momenta J to which individual angular 
momenta are recoupled. Two-particle (non antisymmetrized) states in the two representations are connected through 

|l:(naama);2:(nb/3mb)) =E ^/.™..-.™Jl-('^«");2:K/3); J^^)^ (C22a) 

JM 

\\:{na.a);2:{nbl3)- JM)^ E C'//i^,„„Jl: (n, a m,); 2: (n^ /? mb)) . (C22b) 

nia,mf, 
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The corresponding relations between antisymmetrized states are 

\{na aula); {rib f3mb)) = + S C'/^m„ji,mJ('^a a); (rifc /3); JM) , (C23a) 



JM 



|Ka);(nb/3); JM) = / J] C//4,,„JK am,); ^ /Jm^)) , (C23b) 

V J- + "a/3 ''"ant m^mt 



where the factor -^1 + 5ap Sn^m ensures the correct norniahzation of the antisymmetrized state 
|1: (tt-o a); 2: (tIc a); JM), which is non zero for integer values of J. Antisymmetrized potential matrix elements in 
J-scheme are thus related to those in m-scheme by means of 

yJMj'M' _ \ \ c^JM (-ij'M' Y (C9d\ 

and, conversely, 

Vabcd = Vl + Sap <5ri„)ii, \/l + 5^5 Sn^iia ^ jamajbmb^jam^jdma Kilan6/3Ti<:7"d'5 • (C25) 

JMJ'M' 

Since nuclear potentials are rotationally invariant, they do not depend on M or Ai' and are non-zero only for J = J\ 
such that one can define 

Ki^„a"Jn,7nd<5 = ^J-^' ^MM' yn}ntn^la ' (C26) 

which allows rewriting Eq. (jC25[) according to 

Vabcd = ^ \/l + 5al3 Sn^nt \/l+S^S Sn^na j,,mb C'/^m,jdmd Ki'intnlSd ■ (C27a) 

JM 

Similarly, one has 

^abcd = ^ \/l + 5ai3 Sn^n^ V^l + ^iS <5iic"d Vb^d ^j^ajb-rnb^jam^ja-md ^njntnfna ' (C27b) 
JAI 



JM 



5. Block-diagonal structure of self-energies 

a. First order 

The goal of this subsection is to discuss how the block-diagonal form of the propagators and interaction matrix 
elements reflects in the various self-energy contributions, starting with the first-order normal self-energy 
Substituting Eqs. (|C27aP and (|C19|) into Eq. (|B17p . and introducing the factor 

QnbU.nd ^ ^l + 5ap5^^^, Vl+^75'5«e«. , (C28) 



one obtains 



cd,k 

E\ \ \ ^ f^aricnbrid f~i.JM ^JM •TrJl'^lPl] y^k* -nfifc 
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"cMd 1 J 

^a^-J^^^.AH,^, (C29) 
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where the block-diagonal normal density matrix is introduced through pab = 5ap i^rn„mi, Pnlnt^ such that 

and properties of Clebsch-Gordan coefScients has been used. The fact that the interaction conserves parity and charge 
yields 5^,^^,^^ and Sq^q^, leading to 6ap = Sj^j^ S^^-^^ ^QaQb- Similarly for I]22(i) 



^ab ^ 2-^ ^bcad c d 
cd,k 



^aP'Jmanib J a-jaj 2j -(-1 ^"bncnaWd f^nd?: 



— Uaj3 Umarrii ^n„n(, 

— -S a 8 a'"! 

= -5a/.<5„[AM„J*. (C31) 

Let us consider the anomalous contributions to the first-order self-energy. Substituting Eqs. (|C27bp and ([CT9| into 
Eq. (|B23[) one derives 



^12(1) _ 1 -Vk*fjk 
^ab ~ 9 ^abcd l^c 
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cd,k 
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-A „A i \^ \^ /•no"6nond_- _ / 1 \2.7e V^Jc + ^ t70 [0077! /^H 
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ricTtd 7 V Jo. 

= A « A yi2H(i) 

= 5o.0 <5™„™, /i[rj„, , (C32) 
where the block-diagonal anomalous density matrix is introduced through pab = S^p S^amb Pn}nb, such that 



pi" 



n^nb 



rib m ria [a] "> ' 



It is interesting to note that the first-order anomalous self-energies only involve J = matrix elements as a direct 
result of dealing with a J = many-body state. The other anomalous term is similarly obtained from Eq. (|B24[) and 
reads 



(1) _ 1 T/ Jlk* yk 

^ab ~ 9 Z^ VcdabOt^ 
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cd,k 

\^ \^ \^ fna-ribn^nd f^JM f^JM J [770/?] i/rifc yiifc* 
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n^ndrik 7 mc JA/ 



— "apUm^mb^ J aaj-y "a'^cl, J-j / 2 ■ ^ ^ n^ndU^Ub PudU^ 

TLr^ndUk 7 \ Jo. 

— UafI 'Jniamb ^riaUb 

= SapSm^rabht}i. (C34) 
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b. Second order 



Block-diagonal forms of second-order self-energy contributions ([77|) and ([75)1 can be obtained by considering ex- 
plicitly the angular momentum couplings of the three quasiparticles to their total momentum Jtot, separately for the 
six objects A4, Af, V, Q, TZ and S. One proceeds first coupling particles 1 and 2 to some momentum Jc, which is 
afterwards coupled to particle 3 to give Jtot- The recoupled A4 term is computed as follows 



mim2m3 A/c 
mim2m3A/c rst 
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^ "riantn^n, '^n,. ^n, [K2] ^ nt [1^3] 



'/R3 JQK3K1K2 79 • 11^ JaMajk^rnk^^JaMajk^mk^ 
miManr-n.nt V -^Ja -r 

y T>Jc[aK3KiK2] 7 /"''I 7/"'=2 -i;"'=3 

^ ''nantn^ri, '^n, [ki] '^n, [K2] '^nt [K3] 

— —At • A., \ ^ fnalltrir-ns v'^Jc + 1 ^ -jg f>JefaK3KiK2] 7/"'=l 7/"''2 i;"'''^ 

— ".Jtat3a"Mtotma ^ ^fes^oKaKiKa /2F~rT ' '"^"antrirns '^n, [ki] '^n^ [k2] "^rit [/ta] 

^^jt„t,.'5Mt„t™.A^:r[:::":'.2]j.' (css) 



where general properties of Clebsch-Gordan coefficients have been used. Similarly one derives the M term 



\rk1k2k3 r r \ " --Ttgnf Tigris V^Jc + 1 , -■ x ,7^ + 7^, - ,„ y JJa^aKi K2] i;"^! ■11"''°=^ 7J 

■'^a(J,Jt„t) ~ ^fc3 JQK3K1K2 7 2^- _|_ 2 ^ ' ^nantTlr-n, [ki] [k2] '^i 
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One can show that the same resuh is obtained by recoupUng directly M , as fohows 



\rklk2k3 _ \ ^ f-<JcMc f^JtotMtot \rk1k2k3 



mim2m3_A/c 



— ^ ^ ^ ^ ^ ^ Skip Sm^.^-mr ^K,2(r ^riik^-ms ^K3T Sm^^mt VaVt fa 



I ar pa 

mim2m3Mc rst J^AI^ 



jkirnkj^jk2"ik2 JcMcjk^rnkg ja-majt-nit^jr-mrjs-ms ^ri„ntn^Tia '^rv [p] ''^ris [cr] )it [t] 

'/a '/';3 JaK3KiK2 jki-ni'ki]k2™-k2 JcMajk^rnk^i ja-majk^-mk^ jkirnk^jk2'"i'k2 

mim2ni3Mc n^nsnt J^M^j 

y y-J„[aK3KiK2] y"*:! v;"'=2 //"''a 
^ '^nantn^ris '"^n^ [ki] [/t2] nt [fta] 

\ ^ \ ^ „ „ fnantUr-ris r^JtotMtot /^JcM^ -fyJc[aK.3KiK2] y"-''! W^''2 l/^''^ 

'la. 'Ik3 J aK3K.iK2 ^,JaMajk3mk^^ja-majk3-mk3^rL^ntnrns [k2] rit [K3] 



ES^ n TV, f^o,ntn^n,V^^c^^i^\]^+jk3-Ja(^JtotMtot f^Ja-m^ 
'la " k3 J aK3KiK2 /2T~rT JaMJk3mk^^ JaM^jk3ink3 



V T/'^<:["''3Kl'i2] y"''l y"'=2 Jj"''^ 

^ 'nantn^ris '^n^ [ki] [K2] [«;3] 



A, .A,, \^ „ TT, f»a"t»r»a '^^^g + f 1 V/c -jfc , -.7a T A Jq faK3 «! K2] I j"*" 1 y"l'2 7y"'=3 

"Jtat]a"Mt„t~ma 'la Trk3 J aK3K.iK.2 i r, ■ ^ ^ \ ^) ^n„,ntn^n, '^n, '^n, [sa] "t N] 



which recovers relation (|72ap . The remaining quantities (see Eqs. ((69)) and ((70|) ) are related to M. and A/" by 
permutations of {kx^ki^k^ indices and can be obtained from Eqs. (jC35p and (jC36p by taking into account the 
different recoupling of , jfej and j/cg to Jtot and Jc as follows 



M2k3 ^ y (_i)./c+j.+j.2+j.3 ^2J, + iv/2Jd + 1 ( "^.^Xu^ff j , 

aiJoJtot) Zl^^ ' ^ ^ \jk3 Jtot Jd \ aiJaJtot) 

Jd 



5jt.tjjMt.tnro E E ^'^^ /a«2"«X^ ^^^(2-/<i + 1) (-l)^^+^'=3 | 



jk3 Jtot Jd 



„ T/Jd[aK2KlK;3] 7/"''l 7y"'=3 -i;"'=2 

^ ''n„ntn,n, '^n, [kj] '^n, [K3] '^nt [K2] 



Qk,k2k3 ^ Y^_^y^ + j,+j,^_+j,^ ^ ^ ^ ^ r Jk2 Jfei \ J^l'ik^k^ 

a{JoJtat) Z^^ ' V V 1^ J^^^ J a(,7<iJtot) 



Sjtatja^Mtatm, 



y T/-Jd[aK2KiK3] y"*.-! y"fc3 7y"'^2 

^ "^nantrirns '^rir [ki] "s [13] nt [K2] 

= '5jtoUa'^A/t„tm„ 2„„'[aK3Kl^2] Jc ' (C39) 
a(JcJtct) A^^ ' ^ ^ \]k3 Jtot Jd j a{JaJtat) 

J d 

- -8i ■ 8nr \^ X^TTz f"a"t"r». ^'^Jc + ^ (-9 /-, 4_ -| 'I (_'\\Jd+3k^+3a / j'^i ^^2 -^c 

rLrUsnt Jfi V \. 

V t7J<i[a!'«lK3K2] 7/"''3 7y"'^2 Ij"''! 

^ ''nantn^n, '^n^ [K3] '^n, [K2] "f 

^^j..,jM.„.™aK:^[::;^^::^,^,,^, (C4o) 
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^/Ci/C2/C3 

a {JcJtot) 



Jd 



„ -rA./£i[Q;KlK3K2] -);"'=3 -i;"'»2 7/"'°l 

^ "^nonfrirris '^n^ [ks] Ms [K2] nt [ki] 



™.2 ^27^+1^ '^ ^^ ^ \jk, Jtot Jd 



jki jk2 Jc 

'3 



(C41) 

These terms are finally put together to form the different contributions to second-order self-energies. Let us consider 

^afc^^ ^ example (see Eq. ([75)) ). By inserting Eqs. (jC35|) and (|C36p and summing over all possible total and 

intermediate angular momenta one has 

(\ * ^ 
« rfcife2 ks 



,11(2') 
-"ab 



\ E E 



A^fclfe2fc3 / KAkxk^k^ 
■^^^a(J^Jt„t) l^-'^'h (JcJtot) 



Jtot^totJc k\k2k^ ' 



/3 OmaTnf, 



^E E E 

J "fcinfc2"'=3 ''l'*2'i3 



yM"'=l"'=2"fc3 

rio [aK3K;iK2] J 



■0 {^Zi 



2 "'«3 
[QK3K1K2] Jc 



lW-"'=l"fc2"fc3 

"a [0!K3ftiK2] Jc 



2 ""3 

[aK3KlK2] Jc 



(Wfci +U}k2 +^^3) 



= S ar) yiiH(2') 

— "a/3 'Jmamt, ^n^rit 

Proceeding similarly for the other terms and defining 



rta [ahi3K-ii^2] Jc 



nfcin*:2"''3 

ria [aAi3KiK2] Jc 



1 

1 
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ria [a/«3K,iK2 



Jc ' ra, 



[aKsKiK2] Jc 



^Ua [aft3/ClK2] Jc 



yW-"'=l"'!2"'=3 

ria [QK3K1K2] Jc 



'^Ua [a^3ttlK2] Jc [aH:3ytiK2] Jc 



(C42) 

(C43a) 
(C43b) 



one finally writes 



^l^^'= E E E 



C„ r , , \yni,[o 



infc2"'=: 

"'ria [aK3/Ci/€2] Jc V^^'b [Q!K:3^tlK2] Jc 



rta [cK^t3KiAS^2] Jc/ Tlfa [cK^t3/«iK2] Jc 



(Wfcl -I- Wfe^ -I- LOU 

3) ' 



«^^^= E E E 



n„ [QK3K1K2] Jc \ [QK3K1K2] Jc/ y 



2-j"fci "fc2"fc3 \ ^"fci"fc2"fc3 

Tla [aK3KiK2] Jc y rifo [0^3 i 2 ] Jc 



nfcl"fc2"fc3 Jc K1K2K3 



w - (wfei -I- ujk2 + Wfeg) -I- J?7 



y2lH(2) ^ V V V 



2-j"fcinfc2"*3 /^"fcinfc2"'=3 

ria [aK3KiK2] Jc \ Tlf, [aK3Kin;2] Jc 



w -I- (wfcg -I- Wfei -I- Wfc2 ) - 777 



^nfci"fc2"fc3 ] 25"*! "''a "^3 

ria [aK3KiK2] Jc/ ^fa [aA£3Ki/t2] Jc 

W -I- (Wfc3 + UJki +^k2) - 



5p22H(2) = V V V 

"aMb Z^ 



I? 



ifcinfc2"'=3 /pnfcinfc2"'=3 V Z', 

la [C(K3K1K2] Jc \ n,, [q K3 Ki K2 ] Jc / \ 



^nfcl"fc2"fe3 \ ^"fci"fc2"'!3 

^^a [ati:3ti:iti:2] Jc/ n?, [0^3/^1 /^2] Jc 



nfci"fc2"fc3 -^c K1K2K3 



W - (Wfei -I- U}k2 + Wfca) -I- 7?7 



to + {uJks + Wfci -I- UJk2 ) - 



. (C44d) 



6. Block-diagonal structure of Gorkov's equations 



In the previous subsections it has been proven that all single-particle Green's functions and all self-energy con- 
tributions entering Gorkov's equations display the same block-diagonal structure if the systems is in a 0^ state. 
Defining 



Tab — fJ'^ab — Sap S„ 



rp[a] 



[9a] X 



introducing block-diagonal forms for amplitudes W and Z through 



•\Ajk1k2k3 

'^^fe(Jc./tot) 



JtotifcOA/totTOfc ''^nk [K3K1K2] Jc ' 



ykik2k3 — A., Ti, 57-"''l"'=2"fc3 

^fc (JcJtot) - °Jtot]kOMtot-rnk Vk ^^k [K3K1K2] Jc 



(C45) 

(C46a) 
(C46b) 



35 



with 



'=2 "'=3 ) 7/"fc 



na,a 

and using Eqs. (IC29l) . (fC3T|) . (IC32]) . (IC34l) . (IC44l) . one finally writes Eqs. (gTI) as 



ria [(y.K^n-iK2\ Jc) [ce] 



2] Jc rifc [K3K1K2] Jc 



Ua [q;K;3KiK2] Jc 



*^fcl^fc2^'^3 ^1^2/^3 Jc 

[ l^n^nb M On^nb + ^>-n„nJ l^„, [„] + "riant "n^ [a] 



"fcl"fc2"fc3 K1K2K3 Jc 



QK3K1K2] Jc >ifc [K3K1K2] Jc 



n„ [QK3K1K2] Jc 



life [K3K1K2] Jc 



j.nfcinfc2nfe3 
""rife [k3K;ik;2] Jc 



(C47a) 
(C47b) 



(C48a) 



(C48b) 



The latter four equations constitute the block-diagonal form of Gorkov's equations. Note that pole energies Uk only 
depend on and k, i.e. they display a degeneracy with respect to the magnetic quantum number 771^.. 



Appendix D: ^-functional 
1. Connection between $ and self-energies 

Performing the trace over Gorkov space, the n-th order ^-functional defined in Eq. reads 

*("nG, V^] = - — Ti-n, c ( G" + (") + S^i (») + („) ^22 1 
An I J 

= _^ Tr«, ^ (2 (") G" + (") G^i + S^i (") G^^ \ , 
An I J 

where Eqs. (|46ap and (|49ap have been used to express G^^ and in terms of G^^ and E^^. 
Let us differentiate expression (|Dip with respect to the normal propagator G^^. One finds 



5$(")[G,y] 



1 /"dw' 

~ zln I ^ / 9-77 



JG i(^) 471 1 2^ 



4n 



l^{(4n-2)I]i;(")H + 2S^;("'M} 



(Dl) 



(D2) 



The factor (4r7 — 2) comes from all possible ways of cutting one normal propagation line G^^ in an 7T.-th order self- 
energy diagram Sfifs (") and reconstructing E^^^") by performing the convolution with the multiplied propagator 
G^'^si _ Notice that such a result is not a straightforward generalization of the standard proof of Ref. [7a| as the 
reconstructed self-energy (second line in Eq. (|D2p ') has contributions from all self-energy types (first line in 

Eq. (|D2p '). each containing a different number of G^^ lines. 

Similarly, one can work out the derivative with respect to G^^ , obtaining 



j^(")[G,y] _ 
5G]f{w) " An 



- E 



ab 



doj' 

2V 



^12 {n) 



21 (n). 



SG]f{io) 



Gtlico') + 



SG]f{u^) 



-GllW) 



l-{(2n-l)I]^(")(.) + I]-(")(.)} 



A71 

-^S-(")(.). 



(D3) 
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The factor 2 between the normal and the anomalous case can be intuitively understood as follows. In a closed diagram, 
whenever a G^^ is present, a corresponding G^^ must appear. To a G^^, on the other hand, always corresponds another 
G^^, yielding twice as many possibilities of cutting such a line. 

Summing over n on both sides of Eqs. (|D2p and (|D3|) . together with the analogous ones for G^^ and G^^, one 
recovers Eqs. ()112|) . 



2. Derivation of E^^' from 



At first order, the ^-functional is the sum of two diagrams 



$(i)[G,F]= ioj' 






(D4) 



Using diagrammatic rules outlined in Appendix IB II one can write the corresponding expression 



duj' duj" 



[yacbdG]^^{uj')Gll{u 



ibcd 



duj' 
2^ 



Y.[V,,i,aG]^i{uj')Gll{u')\ . (D5) 



abed 



Applying Eqs. (jll2p and employing Eq. (j46ap one can recover first-order self-energy terms computed in Appendix 
iBTal 



.12 (i; 



(c) 



.22(1) 



= -2 



^Gp:r 

k cd ab 

ab 

cd 



1 f duj' 

2 I 27 



Vabji Gil {^') 



(5$(i)[G,F] 



2SGmu:) 



2<5G/(a.) 



J' 
27 



^"27 

dw' duj" 
^"27 



E[^-..Gii(a.')Gii(a.")] 



abed 



(D6a) 



(D6b) 



(D6c) 



'ibcd 



(D6d) 
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